The order, shape and critical point for the quark-gluon plasma 

phase transition 

Ismail Zakout 1 ' 2 ' 3 , Carsten Greiner 1 and Jiirgen Schaffner-Bielich 1 
1 Institut fur Theoretische Physik, 
2 Frankfurt Institute for Advanced Studies, 
J. W. Goethe Universitdt, 
D- 600 5 4 Frankfurt am Main, Germany 
3 Jefferson Lab, Harvard University, 

Qh' 

<L> 1 Cambridge MA 02318, USA 

OO 1 (Dated: February 6, 2008) 

(N : 

^ ; Abstract 

> ' 

■ The order, shape and critical point for the phase transition between the hadronic matter and 

■ quark-gluon plasma are considered in a thermodynamical consistent approach. The hadronic phase 

m ; 

is taken as Van der Waals gas of all the known hadronic mass spectrum particles mu < 2.0 GeV as 

vO ■ 

well as Hagedorn bubbles which correspond hadronic states with mass spectrum mu > 2.0 GeV. 

43 ■ 

The density of states for Hagedorn bubbles is derived by calculating the microcanonical ensemble 
O ! for a bag of quarks and gluons with specific internal color-flavor symmetry. The mixed-grand 

and microcanonical ensembles are derived for massless and massive flavors. We find Hagedorn 

> : 

bubbles are strongly suppressed in the dilute hadronic matter and they appear just below the 
line of the phase transition. The order of the phase transition depends on Hagedorn bubble's 
internal color-flavor structure and the volume fluctuation as well. On the other hand, the highly 
compressed hadronic matter undergoes a smooth phase transition from the gas of known mass 
spectrum hadrons to another one dominated by Hagedorn bubbles with specific internal color- 
flavor structure before the phase transition to quark-gluon plasma takes place at last. The phase 
transition is found a first order for the intermediate and large chemical potentials. The existence 
of the tri-critical point depends on the modification of the bubble's internal structure specified by 
a phenomenological parameter 7 oc /ig in the medium. 
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I. INTRODUCTION 



The thermodynamical description of the strongly interacting hadronic matter gas can be 
approximated by a free gas with a modified level density. This new level density is given by 
the statistical bootstrap equation and its solution is the asymptotic mass spectral function 
for hadronic fireballs with mass exceeding 2 GeV. A solution to this equation exists only for 
some range of parameters [1]. However, the bootstrap model with internal symmetry [2] of 
the fireball provides subsidiary variables and allows for new types of the phase transition [3- 
5]. The hadronic fireballs are bags of confined quark and gluon components and stand 
for highly excited exotic hadrons; we denote them in the following as Hagedorn bubbles. 
Furthermore, these Hagedorn states have received much attention to understand the phase 
transition in terms of the AdS black hole duality [6, 7] and the gauge field [8, 9]. 

Quarks and gluons are confined due to color confinement while the low-lying hadronic 
mass spectrum is generated by the broken chiral symmetry. The hadron's constituent quarks 
are massive even for light flavors. However, the constituent quark mass decreases slowly with 
respect to temperature and then drops quickly to its current mass at the critical temperature. 
It also decreases with respect to the baryonic chemical potential. The decreasing of u,d- 
quark's constituent mass near the phase transition is not fully understood. Furthermore, the 
internal color symmetry of the bound state remains to be that of a color singlet even for finite 
temperature and chemical potential due to color confinement. When the temperature reaches 
the critical one, color is expected to be liberated and chiral symmetry is restored. QCD 
predicts a phase transition from the hadronic gas phase to a deconfined quark-gluon plasma 
phase. The order of the phase transition for hot and dense hadronic matter, however, remains 
unclear whether it is first, second, higher order, or just a rapid but smooth crossover. The 
internal structure of the hadronic quark and gluon bag remains in a total color singlet state. 
In more realistic calculations, the internal structure is imposed in the partition function for 
Hagedorn bubbles [10]. In the standard MIT bag model, the quarks and gluons are confined 
within a sharp surface, which represent a boundary between different media: a perturbative 
vacuum inside the bag where the quarks and gluon can be approximated as an ideal gas, due 
to the asymptotic freedom of QCD, and a nonperturbative vacuum outside the bag where 
free quarks do not exist due to confinement. 

The hadronic phase consists of the whole hadronic mass spectrum including resonances 
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of all the known particles. Hard core repulsive forces can be represented by an excluded 
volume. Effects from strong interactions are included by adding a free gas of Hagedorn 
bubbles which are bags of constituent quarks and gluons with specific internal color-flavor 
structure. Despite of complexity of the internal color-flavor structure, Hagedorn bubbles 
remain to be in an overall color singlet state. With increasing baryonic density at low 
temperature, the bubble size grows but retains its own internal color symmetry. This means 
that volume fluctuation is expected to be suppressed whenever the bags start to overlap 
with each other for a large chemical potential. When the temperature increases, the surface 
is smeared out until finally the bubbles dissociate at the critical temperature. Therefore, it 
is expected the bubble volume fluctuation increases when the temperature increases. 

The (grand) canonical ensemble and its Laplace transform to the micro-canonical ensem- 
ble for gluonic bags or glueballs was derived by Kapusta [11] without imposing any color 
constraint. The internal symmetry constraint was originally introduced for the statistical 
bootstrap model [2]. The bootstrap density of states can be derived from the MIT bag 
model. In the hadronic phase, the highly excited fireballs derived from the bootstrap equa- 
tion are Hagedorn bubbles (e.g. gas of bags) of confined quarks and gluons in a color singlet 
state or colorless charge. 

The colorless bubble is usually less restricted than the color singlet state. The color 
charges are set to zero for the colorless state uq { = = 0, where Q v is the bubble grand 

potential density and <pQ i are the color parameters in the SU(N C ) representation. In the 
less strict constraint the color charges n 3 and ng are set to zero for the gas of bags instead 
for any individual bag. The asymptotic approximation for the colorless charge is given by 
i(f>Q i — > 0. Hence the colorless bubble is usually approximated to an ideal gas of quarks 
and gluons with the conserved color charge i(f>Q i — * PqJT ~ 0. On the other hand, the 
color singlet state is the ground symmetry projection for the SU (N c ) representation. This 
resemblance between the color singlet and colorless states requires to search for the color 
saddle points around <C 1 and to avoid the Cauchy integral which is used extensively 
in the literature in the large N c limit [12]. 

Gorenstein et. al. [13-15] have studied the gas of bags which correspond Hagedorn states 
using the isobaric partition function. They measured the volume fluctuation for the hadronic 
bags by differentiating the micro-canonical ensemble with respect to the fixed bag volume. 
The micro-canonical ensemble only measures the mass spectral density and does not include 
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the actual volume fluctuation. The real deconfinement phase transition from the hadronic 
phase to the color quark and gluon bags did not take place in their approach. The au- 
thors of Refs. [13-15] argue that the phase transition is only possible when Hagedorn bubble 
expands rapidly to form a big color singlet quark-gluon droplet which occupies the entire 
space. Recently, they have pointed out that the phenomenology of Hagedorn bubble in- 
ternal symmetry decides the order of the phase transition for low chemical potential and 
high temperature [16]. Unfortunately their model fails to predict the phase transition for 
deconfined quark-gluon plasma but instead a new phase of matter appears which retains the 
internal symmetry of Hagedorn bubble. They also were unable to explain how the order of 
the phase transition can be switched from a lower one to higher ones. On the other hand, 
Auberson et. al. [17] have shown that the phase transition to a deconfined QGP is granted 
when asymptotic volume fluctuation is taken into account correctly. They emphasize that 
the micro-canonical ensemble measures only the mass spectrum for a specific bag volume, 
which is not the actual volume fluctuation for the quark and gluon bag. Furthermore, it is 
well known that the bootstrap density of state can be obtained when the bag volume is fixed 
to v — m/AB. Auberson et. al. [17] have approximated the volume fluctuation by allowing 
the bag mass variation with respect to the volume up to the second order in the distribution 
function. Furthermore, they have argued that by relaxing the internal symmetry constraint 
for the color singlet state the resultant Gaussian-like volume fluctuation leads to a second 
order deconfined phase transition. Hence, the color singlet constraint imposed on the bag 
states is not a critical for the phase transition's existence. Therefore, the appearance of the 
deconfined phase transition depends essentially on the bubbles' volume fluctuations beside 
their internal structure constraints. 

Hagedorn bubbles likely appear for the highly compressed hadronic matter. These Hage- 
dorn bubbles coalesce and form hadronic bubbles foam. This state of matter is not a decon- 
fined phase. However, when this highly compressed foam of Hagedorn bubbles is heated, the 
hadronic foam undergoes a smooth phase transition and forms a big quark-gluon droplet as 
the surfaces between bubbles dissolve. At low baryonic density and high temperature, the 
situation is rather different. The formation of Hagedorn bubbles is unlikely to take place in 
the hadronic phase and whenever they appear because of the thermal fluctuations and they 
shall be suppressed by the gas pressure of the external hadronic mass spectrum particles in 
particular those relatively large Hagedorn bubbles. However, when the temperature reaches 
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the critical value, some explosive hadronic bubbles can appear in the system. Whenever 
the internal pressure of the thermal fluctuated bubble approaches the external one for the 
gas of the mass spectrum hadrons, these bubbles grow up and expand rapidly forming the 
deconfined quark- gluon plasma at last. 

The outline of the present paper is as follows. In section II, we derive the canonical 
ensemble for a quark-gluon bubble with internal color-flavor symmetry with massless and 
massive flavors. The micro-canonical ensemble is derived from the inverse Laplace transform 
of the grand-canonical one. In section III, we review the isobaric partition function for the 
gas of bags with excluded volume. We also summarize the conditions of the phase transition 
in the isobaric partition function construction for the gas of bags. The small excluded volume 
for baryons and mesons and large excluded volume for Hagedorn bubbles are introduced in 
section IV. In section V, we present two models for the volume fluctuation in the isobaric 
partition function. The order and shape of the phase transition to quark-gluon plasma is 
analyzed in section VI. Then we present our scenario for the phase transition to quark-gluon 
droplets or plasma for low and intermediate chemical potentials and to foam of Hagedorn 
bubbles for large baryonic chemical potentials in section VII. Finally, we give our conclusions 
in section VIII. 

II. GRAND CANONICAL AND MICROCANONICAL ENSEMBLES FOR A 
QUARK AND GLUON BUBBLE WITH INTERNAL COLOR-FLAVOR SYMME- 
TRY AND MASSIVE CONSTITUENT QUARKS 

In this section, we derive the (mixed-) grand microcanonical ensemble for an ideal gas 
of quarks and gluons which are confined inside a bubble of specific size which carry specific 
internal color-flavor symmetry. Contrary to previous calculations [15, 18-23], we perform the 
derivation also for massive constituent quarks. The microcanonical ensemble is calculated 
by taking the inverse Laplace transform of the mixed grand canonical ensemble. The volume 
fluctuation is introduced by going beyond the standard MIT bag model description with a 
sharp surface to bag models with an extended surface and bubbles with volume fluctuations. 
We assume that the volume fluctuation is measured by the bubble volume distribution 
function. In this context, we propose that the volume distribution function is determined 
by a smeared volume and mass relation relaxing the corresponding constraint of the the 
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standard bag model. 



A. (Mixed-) Grand canonical ensemble 

The grand canonical ensemble for an ideal gas of quarks and gluons confined in a specific 
volume can be calculated by [11, 15, 18-24] 

MP,*) = TrV a e-^™, (1) 

where H is the Hamiltonian of the physical system, N is the conserved quantum number 
operator. The imaginary chemical potential i? can be written in term of fugacity id = log A 
in the mixed grand canonical ensemble where A = e M / T . The operator V a selects those 
configurations that are allowed by specific constraints due to the internal symmetries [2, 15, 
18, 19, 25] of the system 

V a = V r t, =Q -%--- , (2) 

where Vj project the color component such as the color singlet state, Vp transforms the 
system into the center of mass frame [24] while V$ conserves the total number of particles 
in the fireball. In the present model, we consider a quark-gluon gas with specific internal 
symmetry by introducing three projectors. The color projector of the subspace of all states 
that transform under the representation j of the irreducible representation Q(g) of a compact 
Lie group SU(iV c ) reads 



Vj = dj / dfi(g)xj(g), (3) 

JSU(N C ) 

where 

dM = ^ (^) ^ ft ( 2 sin i^T^) ) 2 2n6 ^ 9 n j d9 ld 9 2 ■ ■ ■ d9 Nc , (4) 

is the normalized Haar measure on G(g). However, for small color angles ^<1, the Haar 
measure can be approximated by 

Mg) ~ ^ (^) ' ( II ( e n-e m A 2-k5 (X> t ) de 1 de 2 ---de Nc . (5) 

^ ' \n>m I \n=l I 



It is also convenience to introduce the set 0j where 6i = (pi for i = 1, • • • , N c — 1 and 
4> Nc = — Ylk^^k- We have introduced the set {(fii},i = 1,...,N C — 1 for the group 



representation SU(N C ) while the set {9i},i = 1, . . . , N c is introduced for the group U(N C ). 
In context of set fa the Haar measure reads 

d ^~Jf\\2^) ( II ~ <i>mf I dfadfa-'-d^N^. (6) 

The quantities dj and Xj(flO are the dimension and the orthonormal basis j of the repre- 
sentation, respectively. The basis Xsinglet^) = ^ wr th dj = 1 projects on the color singlet 
state for the quark-gluon bubble. The zero momentum projector for the system enclosed in 
a volume (V) reads 



- _ f <^ iP G -R iP Q -R JPq-R 

" V> (v) e 



f piPPg-f JPPo-r ipPyj-r 



d 3 i 

t(v)/p '' 

where we have replaced R = f3r. This constraint means that we are working in the center 
of mass frame of the quark and gluon bubble [24]. However this projection correlates the 
quarks and gluons in a specific momentum configuration. This projection constrains the 
momenta of constituent quarks and gluon in a specific alignment 

i j k (y) * ' 

The thermal bath breaks the Lorenz invariance. Therefore, it is expected that the high 
temperature breaks badly this projection. Breaking the momentum projection modifies the 
density of states significantly. We retain this projection in our calculations although it is 
not a realistic one for the dense and hot medium as far it is related in somehow to flavor 
correlations. The argument of breaking this symmetry smoothly with temperature is the 
beneath of the introduction of the phenomenological modification of the density of states. 
We relate this modification to the color-flavor correlation. However, in order to soften the 
equation of state for cold quark matter it is also possible to impose a higher order constraint 
in the momentum space such as 

d 3 R G r d 3 R Q r d 3 R z 

< {v) iyT J { v)~(yT J { v) 

On the other hand, for the crystallized quark matter, it is convenience to introduce the 
orthogonal representation O(N) for the rotational symmetries in order to soften the equation 
of state. These aspects is not considered in the present work. 



j k 



The conserved charge, such as the baryon number Nb, is guaranteed by 

% = r^ e «[(iv Q -%)-iv fl ]_ (10) 

J-n 2VT 

The grand canonical ensemble is calculated from the mixed grand canonical ensemble (note 
that as we project on certain quantum states for the bubbles, we are not dealing with the 
grand canonical description in a strict sense; however, the overall dependence of the system 
on the chemical potential is still kept, therefore the term mixed g.c.e.) 

°^e-^Z v {M). (11) 

The Hilbert space of the gas of quarks and gluons has the structure of a tensor product 
of the three Fock spaces for gluons, quarks and anti-quarks. The mixed grand canonical 
ensemble for an ideal gas of quarks and gluons becomes 

Z V (P^) = V [Tr G U G (g)e-^] [Tr^^e^e^] [Tr^^e^e^] . (12) 

The total ensemble is given by the product of ensembles for the constituents gluons, quarks 
and antiquarks. The internal structure for each species is introduced by the representa- 
tion Ui(g). The notation Tr^ represents the traces over the energy states and the internal 
structure representation. Here, we have taken the operator V = V^p_^ ■ Vj, where Vp_ 
selects the zero momentum and Vj projects on specific internal color-flavor symmetry. The 
imaginary chemical potentials can be written in terms of fugacities A« = e %{>i . The analytical 
continuation to the imaginary chemical potential by applying Wick rotation to the conserved 
charge Fourier parameter 

ift = (13) 

introduces the real chemical potentials for the conserved charges in the mixed grand canon- 
ical ensemble. The fugacities for Hagedorn bubble A = \(\ B , \s, ^i) are determined from 
the real chemical potentials Hb,Hs and \xq for baryon, strange and isospin chemical poten- 
tials, respectively, where \b = e^l 1 \ X s = e^ s l T and A/ = e M/ / T . The constituent quarks' 
fugacities are determined by Aq = X 1 ^ 3 and A s = A^A^ 1 for an isospin symmetric fireball. 

The gluons and quarks satisfy the Bose-Einstein and Fermi-Dirac statistics, respectively. 
Furthermore, the gluons are represented by the adjoint color SU (N c ) : N£ — 1 representation 
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while the quarks by the fundamental color SU (N c ) : N c representation. The three traces in 
the Hilbert spaces read [18] 



TrU G (g)e- fi6a = exp 



i-^Tr c ln 

I a 



exp <^ - (V) d G 



d 3 p 



Tr c In 



1 - R adj(^) e 



-Pp 



(14) 



for gluons and 



Tr c U Q (g)e-^e im * = exp <j + Tr c In 1 + R {nnd (g) e -^e + ^ 



]}■ 



ew[+(V)d Q J ^Trdn [l + R iund (g)X Q e-^} } ,(15) 



for quarks, and, finally, 



Tr c %(^)e- /3 %e-^Q = exp I + £ Tr c In 



i + r: 



exp <^ + (V) d Q 



d 3 p 



fund 



Tr c In 



l + Rfond^)^ - ^^ (> 16 ) 



for antiquarks where Eq(p) = \Jp 2 + m Q is the constituent quark kinetic energy and d G = 2 
and dQ = 2 are the degeneracies for the gluon polarization states and quark spin states, 
respectively. The trace Tr c runs over the color index. The matrix R- ac jj(<7) is the adjoint 
color group representation for gluons while the matrices Rf^d^) and R| UIK j(flO are the 
fundamental color group representations for quarks and antiquarks, respectively. 

The flavor index is suppressed in order to simplify the calculations. However, the specific 
color-flavor correlations can be studied straightforward using the present formalism. The 
fundamental representation for the color compact Lie group SU(iV c ), respectively, read 



R fund(^') 



\ 



o e^c- 1 



V 



o o 







( 



) 



,ik<f>i 



Q Q e ik(t>N c -i 



\ 



V 



o o 







, (17) 



/ 



where the set {9i} is suitable for the unitary representation U(N C ) while the set {4>i} is suit- 
able for the special unitary representation SU(N C ). The trace for the adjoint representation 
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reads 



*adj(/)] = Tr c [Rfund(^)] Trc [R*f und (/)] - 1, 



^^cosA;(^-^) + (iV c -l), (i^j). 
i=i j=i 



(18) 



The sum of states is approximated by 

a 

The surface tension 

cPp 



-ci 



d 3 p 
(2^)3 



(19) 



(20) 



is significant for small bubbles and could lead to a bubble instability under specific circum- 
stances while it is negligible for large bubbles. However, the surface effect and curvature 
terms can be absorbed by introducing the so called "quantum volume fluctuation" . We will 
treat the quantum volume fluctuation in detail in Sec. (II C). The mixed grand canonical 
ensemble for a gas of bubbles in the limit of large volumes becomes 

d 3 r 



Z v (P,-&) = dj / dfi(g)xj(g) 

>SU(N C ) J(V)/f33 



exp(Tr c [z Q Q G (r,s)]), (21) 



where 



?qqgM) = Tr c [ z qqgM)] =ZQ(r,g) + ZQ(r,g)+ZG(r,g). 



(22) 



Each integration in the multi-integration space is evaluated by making use of the saddle point 
approximation in the steepest descent method. Using the Taylor expansion, we expand the 
function Tr c [zqqq ( r , 9({@i}))] around the extremum points for the configuration radius r 
and color angles <?o({#i}), 



z qqg(^9) ~ z Q Q G (r ,g ) + 



1 d2z QQG( r >9 



2 dr 2 

N, 



(r - r o y 



ro,90 



1 d 2 z o75a {r,g) 



" QQG\ 

de? 



(23) 
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and 



z qqg(^9) ~ z Q Q G (r ,g ) + - 



1 d2z QQ G ( r >9) 



dr 2 



(r - r ) 5 



+ 1 d2z QQG^9) 



(0i - 0io) 2 ? 



(24) 



»*o,9o 



for the unitary representation U(N C ) and the special unitary representation SU(N C ), respec- 
tively. The Hilbert space for a gas of quarks reads 



Trr 



u Q {g)e 



-f3H Qe iDN Qe i(3p-r 



exp(Tr c [z Q (r,g)]), 



(25) 



where the polar integration over the exponential term is evaluated as follows: 
z Q (r,g) = +2(V) J ^J, 2 In [1 + R fund (»)Age-^C)e'*"] , 

= +2(,/> / ^{-E(- 1 )"s( B *™i(»)v- w » w )"^*" 

- +W / 



oo 



i (Rfmld(j )V— )"^) 

The function S1 °("^) h as an extremum 
(maximum) at r=0 in the configuration 



ra=l 

at r=0 



™ .(26) 
n/3pr j 

unction S1 " r |^^ r ' 1 has an extremum at r=0 and consequently Eq.(26) has an extremum 
imum) at r=0 in the configuration space. In order to evaluate the zero momentum pro- 
• operator, we make a Taylor expansion around the extremum r = for the exponential 
the quadratic term 

(27) 



i "'*Q(3)(r,g) 

z Qm (r,g) » z Qm (r,g) ^ + - — 



1 dh 



where 



* Q (r,9)\ r = = +2(V) J ^ 2 |-2f;(-iri(R fund (^)A Q e-^^) 
= +2{V) J (|^ ln [ 1 + R fund^)- 

The quadratic term 



(28) 



d 2 z. 



dr 2 



r=0 



. reads 



^ j-sf(-l)"i (Rfund(«)^e-" E »«)" 
2(v) [ fP f K V gfundMV^ ' 



(29) 
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71=1 



The above expansions are derived using the following relations 
and 

^{-l) n nx n = -x/(l + . 



(30) 



x) 



(31) 



n=l 



The fundamental matrix R-f unc [(fi , )| 9 o({e i =o}) is approximated to a unit matrix near the ex- 
tremum color angles {9,i = 0}. Hence, Eq.(29) becomes, 



d 2 z Q (r,g) 
dr 2 



= Tr f 



ro,90 



d 2 z Q (r,g) 
dr 2 



n>,9o 



= -2{(V)/^)V Q (m Q (3,X Q ), 



(32) 



where 



(27T) 



[1 + A e-^(f)] 

[e 2 - 



^£ A e r e 2_ m 2.2l3/2 

3 L Q/3 2tt 2 L ql J \ [1 + A Q e- 



(33) 



Hence for quark and antiquark it becomes 

de 



r 2 2 n2~\ 3 / 2 

L e " m q P \ 



[1 + Age-f [l + A^e-]" 



.(34) 



For zero quark mass limit, it reduces to 

N r . /7tt 2 



2tt 2 



In 2 A, 



(35) 



On the other hand, the color state projection is calculated by evaluating the integral over 
the Haar measure. The quark exponential term Zq (r , g) around the point r reads 

Z Q (r ,g) = e W (Tr c [z Q (r ,g)])=Tr c U Q (g)e-^e im Q, 

= exp |+2 ((V)/P 3 ) jT ^ 2 e^je 2 - m 2 Q [3 2 Tr c In [l + R innd (g)\ Q e^} J , 

= exp i +2 ((^)//3 3 ) 



de l (2 m 2o2N3/2nn f R fund(^) A O e 
? -^o( e ~ m QP J ir C 



[ 1 + R fund(^)V" e ] 



(36) 
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The exponential term for system of quarks and antiquark becomes 



Z QQ( r o,g) = exp [z Q Q(r ,g)] 



where 



The color internal symmetry for the quark and antiquark is determined by 

R fund(3)V" £ R fund^ A Q le ~ £ 



J Qo( e . 9) = Tr c 



[ 1 + R fund(^) A Q e ^] 



1 + R fund^ V 



and the straightforward calculation for the trace leads to 



ReJ QQ( e >^) = 

i 

= E 



1 + A e^e~ e 1 + An 1 e _i9 < 



Aq 1 c £ cos^j + 1 



+ 



\r>e e cos9i + 1 



A Q 2 e 2e + 2AqV cos ^ + 1 \ 2 Q e 2t + 2A Q e e cos 0< + 1 
The variation with respect to 6*j in the C/ (iV c ) representation is given by 



A Q e- e (1 - Afo-*) 
[1 + A^e~ 2e + 2A Q e- £ cos0j : 



sin 



sin0». 



[1 + A Q V 2e + 2AgV e cos ^] ' 
In the limit & (^-) <C 1 and cos#j 1, we have also for SU(N C ) representation 



d . , 



Age- (l - Aje"*) + A Q V £ (l - A Q V-) 



[1 + A^e 



2 c -2e 



2A Q e-f [1 + A^V* + 2A Q 2 e-] 



x 



N c -1 



sin(0i) + sin( ^ <^ 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



On the other hand, the Taylor expansion of the gluonic exponential part around the 
extremum (r = 0) reads 



z G (r,g) w z G (r,g)\ r=0 + - 



I d 2 z G (r,g) 
2 



dr 2 



(43) 



r=0 



13 



The first term is given by 

z G (r,g)\ r=0 = -2(V) 

= +2<y) 



d 3 p r 

j^y 3 ln [ 1 - R adj(^) f 

d 3 p Pp R adj(#) e 



-Pp 



-Pp 



(2tt) 3 3 



1 - 



R adj(^) f 



-Pp 



The second term reads 
d 2 z G (r,g) 



dr 2 



-2(V) 



r=0 



d 3 p ( (3 2 p 



(2k) 



R adj^)< 



-Pp 



1 -R 



adj(s> 



-Pp 



The straightforward calculations of the trace gives 



z G (r ,g) = Tr c z G (r,g)\ r=0 , 

d 3 p Pp 



+2{V) 
+2{V) 



(2tt) 3 3 

d 3 p Pp 

(2tt) 3 T 



Tr f 



E R adj(sV 



k(3p 



_fc=l 

oo / N c 



5]cosA;(^-^) + (7V c -l))e- 

k=l 



k/3p 



(V) fp 3 d 3 ppp e~^{ 



COS 



6j) - e 



-Pp\ 



P 3 J (2tt) 3 3 [1 - 2e-fo coa(6i - 9j) + e- 2 to] 



+(N C -1)2- 



(V) fp 3 d 3 ppp 



-Pp 



p 3 J (2tt) 3 3 l-e-Pp' 



The first derivate reads 



dz G (r ,g) 



= +£/l(*n-0;), 

inl- 



and 



dz G (r ,g) 
d(f) n 



(Nc-1) 



for U (N c ) and S77 (iV c ) representations, respectively, where 

t (V) f P 3 d 3 p e-fosmil) 



AW = "2- 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



p 3 J (2tt) 3 [1 - 2e~Pp cos(^) + e" 2 *] ' 
The saddle points for the multi-integration over the color parameters are found by maximiz 
ing the exponential term for the quarks and gluons, 

d 

[z G (r ,g) + z Q Q(r ,g)] = 0, i = 1, 



,N C 



d9i 
d 

[zg(to, g) + z Q Q(r , g)] = 0, i = 1, • • • , N c - 1 



(50) 
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for U(N C ) and SU(N C ) representation, respectively. The solution of this set of equations 
is go (0i • • • Nc = 0) (g (0i • • • <f) Nc -i = 0)). The sets (9 U ■ ■ ■ 9 Nc = 0) and (0 1; • • • 0^ = 0) 
are solutions for the saddle points in the U (N c ) and SU (N c ) representations. These solutions 
are satisfactory under the assumption that the angles — f— ) <C 1 are the dominated range 
for the saddle points in the multi-integration over the color parameters. However, in general 
case this solution gives a good approximation for the integration over the Haar measure 
but this approximation does not necessary work correctly for quark-gluon bubbles in the 
extreme conditions in the entire T — fi B plane. Under certain conditions in particular in 
very high temperatures, one or more of the color parameters becomes — 7^ in particular 
at very high temperatures. When such a case takes place, it is expected that another class 
of solution emerges and modifies the density of states. 
In the limit g = go((fii — 0, i — 1, • • • , N c — 1), we have 



+ 



1 + \ Q e" 1 + A, 



and 



1 d 2 
Re 2 W, Xq ^ €,9) 



90 



Q2 

90 

1 - \ Q e~ 



def QQ{ 

Xoe~ e 



A; 



[1 + A Q e- 



+ 



(1 



We define the following terms 
A Q Q (m Q [3, X Q ) 



poo 

J rriQfi 



^llf e 2 _ m 2 / o2x3/2 

2tt 2 3 1 qP > 



[I + A^e-] 



T QQ^9)\ gQ 



and 



poo 
Jm Q P 



de 1 



2tt 2 3 



dffi 



2^QQ^ t ^9, 



no 



> 0. 



In the limit of massless quarks the above terms reduce to 

7tt 2 



A QQ ( m Q = Mq) 
B QQ ( m Q = 0,A Q ) 



Nc 
3 
1 

6^2 



In 2 A, 



120 4 
(31n 2 A Q + 7 r 2 ). 



1 + 



I n 2 A, 
2tt 2 



(51) 



(52) 



(53) 



(54) 



(55) 
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At the color saddle points g ({0* = 0}), the gluonic terms reduce to 

XV) f(3 3 p 2 dp(3p 1 



z G (r ,g ) = +[N C (N C - 1)]2^/ J 
+ (iV c -l)2 



2tt 2 3 [ePP-1] 
{V} f(3 3 p 2 dp(3p 1 



f3 3 J 2n 2 3 [e^-1]' 



2 {V) A 

2 ~p G ' 



where 

The second derivative of the gluonic part with respect to r in the limit #o reads 



d 2 z G (r, g ) 



r=0 



,<!>. 

/9 s 



<9r 2 
where 

7T 2 

£> G = 4(iV 2 - 1) — . 
Furthermore, the second derivative for the gluonic part reads 



and 



dd m dd n 



— ' ' " ' ' - + ^/ 2 (0n - <f>j)6 nm ~ h(4>n ~ <frm)\ m ^n 



d(pmd(j) n 

3+n 



N c -1 

+ /2(-0JV c + 0j) + /2(0n - 0iV c ) + ^(^m ~ 07V C ), 



for U(N C ) and SU(N C ), respectively, where 



= (V> [pWp [(l + e- 2 ^)cos^-2e-^] 
J21W " /? 3 7 (2tt)3 [l-2e-*cos(^) + e- 2 *] 2 ' 

In the limit of the color saddle points g (fa — 0, % — 1, • • • , iV c — 1) we have 

f (3 3 p 2 dp e~P p 



MO) = -2^ / 
= _ 2 (V) # G 



2tt 2 [l - e-Ppf 



P 3 2iV c ' 
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where 



B, 



G 



(64) 



The resulting quadratic expansion for quarks and antiquarks reads 

(V) r 

z Q -Q(r,g) « + 2— [%(mg^ A Q ) 



1 ^ c 1 



i=l 



(65) 



The same procedure can be done for gluons. The expansion for the gluon term around the 
extrema (r , go) reads 



z G (r, g) » 2 ^ - ^ G £ (9? - lp G r 



(66) 



where 



The mixed grand-canonical ensemble is calculated via 



(67) 



r r (p r 

Z V (P, 9) = by ((3, A) = dj / dfi(g)x(g) / 77777^ ex P ( Tr c [ z qqg]) > ( 68 ) 

J5C/(Ar c ) -/<V}//3 3 \ V IIP 



where 



-QQG 



Q 

= 2{{V)/f) [A G + J2(^ Q Q( m Q^ X Q)) 

. \ Q 

- \ \ v g + Yl V QQM Ag) V 2 - i ( + E %?(™Q^ A Q ) ) £ 61? 



.(69) 



Q / \ Q / i=l 

Each integration of this function is evaluated using the saddle point approximation. The 
color singlet state for gluons and flavorless quarks blob is integrated over the Haar measure 
with the help of the following relation, 



c 



ISU(Nc) 



d/iQ?)x s inglet(^ = ex P _ o" ^2° 



r<(i-N 2 )/2 JVc-l 

II A (70) 



.1 = 1 



1=1 



(2vr)( iV - 1 )/ 2 ViVe 



j'=i 
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where 



C 



Bo + E (*< 

Q 



[m 



(71) 



In order to fit the phenomenological prediction about the existence of the critical point, 
we shall generalize this relation for N c = 3 to arbitrary color symmetry configurations with 
the phenomenological parameter 7, 



L iN / Mxs ^ {g)s {P : = °) exp HI?* 2 ) k ^ 



00 



'4(7-1/2) 



,(72) 



where 7 = 1/2 and 7 = 3/2 for colored and color singlet states respectively The above 
result is motivated by using the general result for Gaussian integration 



J dOe^e- ce2 ~C-^ +1 )/ 2 ,Coc 



00 

03 



and 



(1)' 


Id*'™* " 


2 









^det^M 

2 aet de t de k 



-cf(d ) 



(73) 



(74) 



where i,j = 1, • • • , JV C — 1 and /, are the exponents of the pre-exponential. The param- 
eter 7 depends on the underlying overall symmetry of the system and increases when the 
symmetry configuration becomes more complicated. Furthermore, we assume that this phe- 
nomenological parameter 7 also absorbs the breaking of the Lorenz symmetry invariance of 
the momentum projection due to the thermal excitations 

(V) 



I 



drr^e 



(75) 



The momentum alignment, Pq + Pq + Pq = 0, which is introduced by Kapusta[24] in 
order to reproduce the bootstrap density of states is no longer preserved in the extreme 
temperature and diluted matter. The phenomenological parameter 7 increases as the color- 
flavor correlation of the system increases and the color configuration symmetry becomes 
more intricate. We will use this parameter 7 to study the impact of these color-flavor 
correlation on the phase transition to quark-gluon plasma. 
The zero momentum projection is calculated using 



/ 



d 3 r 



■ exp 



^2 D = 



2vr\ 3/2 

(v)/p \d) 



(76) 
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where 



D* = 2 (v G + V Q Q{m Q (3, \ Q ) j . (77) 



The resulting grand-canonical ensemble becomes 

-[4 7 +l/2] 



where 



Z (A (V)) ~ ^ ( ^ ) -3/ 2(C , ) -4b-i/2] (W^ ~ 7 exp (|^cr(^)) , (78) 

d*{(3) = 3 • 2 ^4o + 5>qq(™q/?, A Q ) j . (79) 
B. Microcanonical ensemble 

The asymptotic density of states for the microcanonical ensemble for large fireball energy 
W is calculated by taking the inverse Laplace transform for the grand canonical ensemble. 
The inverse Laplace transform for large W — > oo is evaluated using the steepest descent 
method, 

Z(W,(V)) = — / d(3e^ w Z((3,(V)) 



/(/? S )e (V0/(/3 S. (80) 



27ry (v)f»(/3 s y 

where /3g is the saddle point. The pre-exponential function in Eq.(78) reads, 



gtf) = H^( D *)-S/2 ( ^ ) ^7-l/^ (<V)) -[W/2| ) g3[*y+l/2| j (gl) 

V3 

while the exponential function reads, 

/(£) = ^ + 3^3^^) = f 3x+ 3^^^)" ^ 82 ) 
In order to separate the volume fluctuation from the mass fluctuation, we define the energy 
density parameter x = W/ (V). The second derivative reads 

d 2 f{(3) 



dp 



2 ' 



= 4^-2^ + -^ (83) 
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Eq.(83) is evaluated with the help of Eq.(79). The gluon terms read 

d 2 d 

W G = W G = 



(84) 



The quark and antiquark terms read 
d , , „ , , „ f°° de 



8(3 



poo 

AQQ(m Q P,\ Q ) = -mp I 



2tt 2 



5 2_ ^2)1/2 



[1 + A Q e-] [1 + AqV'] 



(85) 



and 



d 2 f 



de (e 2 -2mp 2 ) 
mQ p 2tt 2 (e 2 -m 2 q P 2 y/ 2 



Age 



A^e" e 

+ 



[1 + A e- e ] [1 + A^e-] 



(86) 



The value of the saddle point /3g given in Eq.(80) is determined by maximizing the 
function f(/3), 



f'(P)\ 



0. 



(87) 



This point has a nonlinear solution for massive flavors 





1/4 




1/4 


X 




3 d*((3) _ 





(88) 



0c 



For massless flavors, the saddle point has a trivial linear solution [15, 18]. Finite constituent 
quark masses do not seem to effect the order of the phase transition although it modifies the 
nonlinear solution for the micro-canonical ensemble. The nonlinear solution of Eq.(88) is 
obtained numerically. The micro-canonical ensemble with massive flavors will be studied in 
our forthcoming work in particular to study quark and gluon bubbles with strangeness. The 
incorporation of strangeness in the canonical ensemble is particularly important for high 
temperatures [26]. 

The remainder of the article will focus on quark and gluon bubbles (i.e. fireballs and 
hereinafter we abbreviate them as FB) with massless quarks. For massless flavors, the 
function d*((3) = d*(0) becomes independent on (3 and the saddle point solution for the 
steepest descent method is known analytically. In this case the inverse Laplace transform 
becomes 

' /3q— ioo 



i pfSo+ioo 

— / dpp n e v ^ x+d ^ 0)/3 ^ 



1 



1 



2^V l l 2 x 1 ' 2 



[d*(0)/x] (1+2n)/8 e^K(°)* 



,1/4 



(89) 
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The density of states a(W, V) aS y m w Z(W, V) with respect to the bag energy W with 
massless flavors reads 



^(WWasyrn = AV^x^- 1 exp [rf* (O)* 3 ] 
where x = (m — BV)/V and 



31 i/4 



« = cf (0), 



3-2 



(iv c 2 - 1 



7T 



90 



v^iV c /77r 2 ln 2 A r 



Q 



120 



1 + 



2tt 2 



The pre-exponential coefficient reads, 



v3 

= Vt 4(7 ~ 1/2) [^(o)] 37 -\ 

8 



where 



(90) 



(91) 



(92) 



D* = -cT(O), 



= 2- 



3 



(93) 



We have introduced the variable x in order to simplify the equation when we study the 
isobaric partition function below. In the MIT-bag like models, the energy for the quark and 
gluon bubble W is fixed by W = m — BV . However, the ansatz for color-flavor correlations 
alters the quark and gluon bubble internal symmetry and have a strong impact on the phase 
transition diagram, as will show below. 



C. Bubble volume fluctuation 

The quantity cr(P, (V)) with cr(P, (V))d A P describes the density of states of noninteracting 
quarks, anti-quark and gluons confined in a bag with the total 4-momentum in the interval 
(P,P + dP) and with the specific volume (V). The canonical ensemble measures only the 
mass spectral density. Therefore, the density of states given by Eq.(90) does not measure 
the bubble volume fluctuation. In order to measure the bubble volume fluctuation we need 
to find the wavefunction for the quarks and gluons blob bound state. The density of single 
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particles levels can be calculated using the multiple reflection expansion [27-29]. The volume 
fluctuation is considered by smoothing the sharp boundary of the MIT bag due to extended 
potential profile. The density of states can then be calculated for a bag with an extended 
surface profile. 

Let the system comprises a gas of bubbles with various species. Each species type is 
labeled by the index [I] and is classified by the bubble internal structure of constituent 
quarks and gluons. The square bracket is to distinguish the species type index from the 
quantum numbers of the bubble's constituent particles. The eigenstate for a constituent 
particle i in the bubble of species type [I] is given by the energy eigenstate E^ ni and the 
wavefunction ^/^ n .(Xi). The label rij indicated the energy quantum number which can be 
occupied by the particle %. In the standard MIT bag model, the interaction between the 
particles is neglected and the confined constituent particle wavefunction for a bag with a 
sharp boundary is given by ty[q ni (X) = N^nji (t^ 7 ") ^ (r) where the Bessel function^ 
vanishes at the bag's surface r = R^, i.e. ji (x ni i) = 0. In this example the bubble's species 
type is specified by the size of the bag radius ify]. 

The level density of states for a bubble of species [I] with a sharp surface boundary can 
be calculated from 

States 

a {I] (E) = V { i] ni 5 {E - E [I]ni ) , (94) 

rii 

where T>^ n . is the constituent particles degeneracy factor of the energy level n« while E^ n . 
is the energy eigenstate and [I] is the bubble's species quantum number. The sum over rii 
runs over all the energy states. However, the level density of states for the bag with an 
extended surface and properly adjusted potential reads 

P[I] (E,r) = J2v [I]ni S (E-E [I]ni ) |^/ K (r)| 2 , (95) 

rii 

where ^[r\m(f) is the bubble's constituent particle wave-function. After integrating the 
configuration space distribution function, the density of states becomes 

<J[i](E) = J d 3 r P[I] (E,r). (96) 

When the system consists bubbles with several species, the density of states becomes 

Types 

p(E,r) = PV\( E > r )> ( 97 ) 

[i] 
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where each species is labeled by type [I]. The (mixed-) grand canonical ensemble for con- 
stituent particles in the bubble of type [I] with specific internal color-flavor structure, reads, 



Xre -/3%, =Tr([/]|e-^|[/]> 
^2 {[I]\ni,n2---) (^i,«2---|e _/3 ^|ni,n 2 ---) (n u n 2 ■ ■ • 



ni,n 2 - 



= |tof (| (I)* |fl (I)II^e-^. = |rf 3 X|M/ [7] (X)| 2 nE e ^ lIK ' ( 98 ) 

i n; " i Hi 

where the sum runs over the occupation number rii for one-particle state while the product 
Y\j is the number of states which can fill the one particle state. The coordinate X is the 
cluster's configuration space while 



constituents 
([/]K,n 2 ...> = *[/](*) I] ^' 



(99) 



is the cluster wavefunction. The function ^f^(X) is bubble wavefunction in the configuration 
space while ni are the constituent particle's Fock states. 
The number of "one particle" states reads 



states 



=1 rii 



[J] 



cxp 



= exp 



states 

ln^e"^ 



x 



d 3 p 



(100) 



The resulting ensemble for quark and gluon bubble of type [I] with the extended surface 
becomes 



Hi] 



|^ [/]( X)| 2 exp {v x J^ln^ 



-P E [I\n 



(101) 



where V x = ^-irX 3 . The density of states for single-particle levels has been studied in [27- 
29] extensively using the multiple reflection expansion method. When we go beyond the 
standard bag model with a sharp surface to a smooth one, it is possible to express the bag's 
volume fluctuation as follows 



SV = Vx-(V m ), 
23 



(102) 



for bubbles with a conserved number density. Therefore, the volume distribution function 
for bubbles with smooth boundaries can be found by the following replacement 

/ «^]» = / dV x P {I] (V x - (V {I] ))f(V x ) , (103) 

which reduces to 

P[/](^-<^]» « 8 {V x - (V {1] )) , (104) 

for bubbles with sharp surfaces. Hence the density of states for a bag with a smoothed 
boundary can be written as follows 

W 

p(W, v) = P (x, v) a(x, w)asym, and x= — , (105) 

v 

where P[/j (V x — (Vj/j)) — > P (x, v) and V x — > v. For the standard MIT bag model, we have 
(V {1] ) = i| and 

P[i] (V x - (V {1] )) ~5{V X - m/AB) = 5 - SB^j , (106) 

where m and W are the mass and energy, respectively, while V x and B are the volume 
and bag constant, respectively for the specified bag. In order to fit the phenomenology, 
the spectral density is the average density for quark and gluon bubbles with smoothed 
boundaries and their volumes are close to (Vj/j). This density of states is given by 

a (e, (V {1] )) = J dV x f Av (V x - (V {1] )) a(e, V x ), (107) 

where the function f& v is a smoothing function normalized to 

' dV x f Av (V x -(V {I] )) = 1, (108) 



and e = E/V x is the bubble's energy density. In the phenomenological calculations, the 
Gaussian smoothing function is given by [27-29] 

"A 



6(x)^f Av (x) = \l^e-^ 2 L« n 



E* 2 



(109) 



The L°[x] denotes the associated Laguerre polynomials, A is the smoothing parameter, and 
E is the bubble's energy (or mass) while (Vj/j) oc E in MIT bag model. The smoothing 
function is analogous to the volume distribution function defined in Eq.(106), 

P[fl {Vx ~ (V {1] )) = f Av {V X - (V {1] )) , 

= fA v (V-V ). (110) 
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In general, the cluster's wave-function is very complicated and in principle is not known in 
detail. However, the standard MIT bag model can be extended by smoothing the sharp 
boundary condition. The bubble with an extended surface boundary is approximated by 
replacing the delta function with a smoothing function as given in Eq.(109), 

S(V - V ) -> f Av iV - V ) . (Ill) 

This can be done e.g. using the multiple reflection expansion [27-29] for calculating the 
density of states and the Strutinsky smoothing method [30, 31]. The volume fluctuation to 
lowest order for the Gaussian smoothing function arrives then at 



5(V-V ) -> \/-*— 



A / E „-%(V-Vof .l A / E 



IT V 7T 



e -#(v-Vo) 2 ? ... (112) 



for n — 0, 1 and • • • , respectively. Nevertheless, the order of the volume fluctuation plays a 
crucial role in the order of phase transition at the extreme conditions. The parameter A is 
chosen small as possible, but nevertheless sufficient large to taken into account the extended 
surface of the bag effectively. The high excitation of the smoothing function f& v (x) due to 
the thermal excitation presumably leads to the bag instability. 

Note that our model differs from the Strutinsky smoothed density [30, 31]. We soften the 
density of states by allowing the volume fluctuation with a specific energy density instead 
of smoothing the energy spectral function with respect to the bag energy. Our ansatz 
can be justified by the multiple reflection method where the potential is deformed by 5r 
and the particle wavelength is varied by L = Lq + SL. A more detailed analysis will be 
considered in forthcoming work to study the volume fluctuation due to potential deformation 
and smoothing the bag boundary by considering a bag with an extended surface using the 
multiple reflection method [27-29]. The Strutinsky smoothed density of states a(E) is 
obtained from the density of states, a(E), by the Gaussian smoothing 

a(E) = J dE'f AE (E' - E)a(E'), (113) 

where f&(x) is given by Eq.(109) but instead Ay = A/E we have A E = A/V . However, 
it is possible to derive the smoothed density for the bubble's volume from the Strutinsky 
smoothed density by substituting E — eV and dE = edV . 
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III. THE GAS OF BAGS WITH VAN DER WAALS REPULSION 



The partition function for N relativistic noninteracting particles reads [14, 32] 

oo oo 

Z(T,V,X B ,Xs) = X b B E *s S Z(T,V,n B ,n s ) 

ng=—oo ng- 
oo _. „ N 

Af=0 J i=l 

N , N 

A 1 J A 1 



3=1 



3=1 " 3=1 



where 



r(mi,vy dw ., • • •) = r(m;,T>y dw ., A B , A 5 , •••) 

oo 

= AnSr ( m ^ T > u VdWi' nB ' n 5'"') ( 115 ) 



riB=— oo 



is the particles density of states with the fugacities \ B = e~ and \s = e t and n B ,ris 
are baryonic and strangeness densities and are the Fourier modes with i6 B = fi B /T and 
ids = fis/T, respectively. The density of states given by Eq.(115) consists the known 
hadronic mass spectrum particles and the continuous density of states for Hagedorn bubbles 

Baryons 

T(m i ,T,vy ( ptf.,\B,\s,---)= Y D F j) i (m,T,v,\i)8(m-m i )5(v-vy$w H ) 

i 

Mesons 

+ Y ^BEi ( m > v i — rrii)5(v — ^vdW//) 

i 

Types 

+ J2 $>/K7>,A/), (116) 

7 rre,i> 

where the sum runs over the baryons and mesons mass spectrum and they satisfy 
Fermi-Dirac and Bose- Einstein statistics respectively. The terms -Dpj^ (m, T, t>, A«) and 
-^BEi ( m ' t '' ^») are ^ ne degeneracies for the Fermi and Bose particles, respectively. The 
continuous densities of states run over the various species of Hagedorn bubbles and for large 
masses they obey Maxwell-Boltzeman statistics. Each species is labeled by the type [/] 
and characterized by specific quantum numbers such as mesonic fireball, baryonic fireball, 
• • • , etc. The excluded volume is taken as the hard core Van der Waals repulsion volume 
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v VdWi = ^ or a s y s ^ em °f identical particles with size V{. The grand partition function 
for a gas of identical particles of mass mi in a volume V with specific quantum statistics 
reads 



Z Stats( T > v \ m ^ A *) = ex P [^Vstats^ m *> ' 



(117) 



where rrii and Xi are hadronic particle mass and fugacity, respectively, while T>i is the degen- 
eracy factor stemming from the internal degrees of freedom. The volume V is the hadronic 
particle's Van der Waals excluded volume. The subscripts Stats = MB,FD or BE corre- 
spond to Maxwell-Boltzmann, Fermi-Dirac and Bose-Einstein statistics, respectively. The 
function ^Stats^' m »' A ») rea ds 



and 



d 3 k 
(2^)3 

d 3 k 



In 



l + A e -v^R/T 



VBE (T; m z , Xi) = -ViJj^ ^ [l - A* V* 5 ^] 



After integration by parts, they become 



^Stats( T ^-^) = ^ deaL 



~T>i I dk 



^Jk 2 + rri; 



where the quantum statistics distribution functions read 



(118) 



(119) 



(120) 



: /Stats( T ; m ^ A ^)' ( 121 ) 



122) 



and 



fpi)(T; ,77ii, Ai) 



/BE( T ! m " A *) 



A^exp^F + mf/Tj+l] 



A^exp/^ + mJ/T) -1 



(123) 



(124) 



for Maxwell-Boltzmann, Fermi-Dirac and Bose-Einstein statistics, respectively. The pressure 
and energy density for each particle are determined, respectively, by 



Pi = T ^stats( T ; m ^ A 0> 
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(125) 



ei = T #r^stats( T ; m ^ A *)> ( 126 ) 

The baryonic density for each particle reads 

d\ d 

n Bi {T-m h \) = T^-—^ tdXs {T-m i ,\), 
91n\ f d 3 k 

/d 3 k 
^/Stats( T ^-^)- (127) 

The above quantities can also be calculated for the anti-particles be simply replacing Aj — > 
A" 1 . The grand canonical ensemble for a gas of non- interacting multi-particle species obeys 
the relation 



Z(T, V- mi , A ls • • • , m n , A n ) = J] Z(T, V; m„ A,). (128) 

i=i 

The isobaric partition function is calculated by taking the Laplace transformation of the 
grand partition function [17, 32, 33] 

dV exp (sV) Z(T, V; m { , vy dWn , A ») 

VdW 

= l/[s- ¥ » Sj (T;m i ,\*)], (129) 

where ^vdW/f * s ^ ne hadron's Van der Waals excluded volume. The effective fugacity 
becomes 

K = exp (-vy dWH s) \ u 

X = exp {-vy dWH s) A" 1 , (130) 

for quark and antiquark, respectively. In order to simplify the notation for different statistics 
cases, we use the "Maxwell-Boltzmann"-like notation 

exp (-^VdWi s ) ^ (T; m u Xi) -> <p s . (T; m h A*) 
^exp (-^vdWi s ) VsA T 'i m h Ai) -> ^^(T; m 8 , A*). (131) 

We have assumed the Van der Waals excluded volume ^vdW// * s ^ ne same for all the known 
mass spectrum hadrons. The isobaric partition function with the isobaric ensemble (T, s) of 
a system is characterized by the external pressure p = Ts rather than the fixed volume V. 
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The point s = s* is the extreme right singularity point in the limit of infinite external volume 
V — > oo. The isobaric partition function with multi-particle species reads [13, 14, 17, 32] 



- ^exp (-*VdWi s ) ^(T; m h A*) 



,(132) 



Z(T, s; m 1 ,vy dWv ■ ■ ■ , m n , ^vdwJ = V 
where 

exp (-^vdWi s ) ^(T; m h \) = f H (T, A; s) + f Q (T, A; s). (133) 



t< 



The first term fii(T, A; s) denotes the ensemble for the known hadronic mass spectrum 
including their resonances and antiparticles in a hot and dense medium, 

n 

f H (T,X;s) = J2[psAT;m l ,Xt)+ VSt (T;m t X)] , (fermions), (bosons). (134) 

The masses of these particles are taken as listed in the particle data group book [34]. In our 
numerical calculations, we have included the spectrum of all 76 nonstrange mesons and all 64 
nonstrange baryons and their antiparticles as well as Hagedorn bubbles for the highly excited 
hadronic states. The effect of strangeness will be considered elsewhere. These Hagedorn 
bubbles naturally appear as fireballs in the heavy ions collision at high temperature. The 
second term in Eq.(133) corresponds to the ensemble of Hagedorn bubbles which exist in a 
color singlet state [14, 15, 23, 32]. The isobaric ensemble for Hagedorn bubbles becomes 

"OO r-oo 



/»00 /"OO 

f Q (T,X;s) = dvdme-^ s p(T,m,v,X Q )if Q (T;m) 

JVq J mo 

oo poo 

dvdme~ 4vs p(m,v)ipQ(T;m), (Maxwell-Boltzmann), (135) 



' Vo Jm 

where p(m, v) measures the mass spectral density and the volume fluctuation. The asymp- 
totic behavior m ^> T of the function <£q(T] m) reads 

d 3 k 



V? Q (T;m) = J ^-^exp (-y/k 2 + m 2 /T^j , 



^] K, {m ,T) 

mT\ 3/2 

f ) -'" /T <*> 

Furthermore, in addition to the bubble's quantum ground state, it is possible to take into 
account the higher quantum excitations of Hagedorn bubbles such as states with higher 
angular quantum momenta as well as mesonic and baryonic dominated bubbles (fermionic, 
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bosonic fireballs), . . . etc. The density of states for the standard MIT bag model (with a 
sharp surface) reads 

p(m,v) = 5(m-ABv)a(W,v) asym , (137) 

where the gas of quarks and gluons is projected on a color singlet state and is confined 
in a spherical cavity with a specific volume v = m/AB. The subscript asym denotes the 
asymptotic density of states for the excited hadronic states. 

The isobaric ensemble Eq.(135) for Hagedorn bubbles with sharp surfaces and density of 
states given by Eq.(137) and Eq.(90) can then be written as follows 

3/2 



e -m,T 



J 

Jv 



Vo 

oo 



roo j-oo /mT\' 

f Q (T,X:s) = J v 5(m-ABv)a(W,v) eiSym ^—j 
-(3^)asym {^j'' 

POD 

= C dvv-^-* /2) - l e- A < s - S0 \ (138) 

JVo 

where W = m — Bv and the prefactor 

C = A{W) 01 - 1 ( ^-M . (139) 

Note, that we have introduced here the parameter 7 which describes the internal color-flavor 
configuration of the bubble (see section IIA). The quark and gluon bubble's internal isobaric 
pressure is given by 

The pressure stems from the isobaric pressure at P = Ts. The external pressure p Q = 
T/q(T, A : s) for Hagedorn bubbles is written as follows 

/ X 47-3/2 -oo 

/q(T,A: S ) = C^|J J z dzz-^-W-'e-' (141) 

where 

z = 4V (s - s ) • (142) 
In order to analyze the phase transition, it is useful to introduce the function 

/oo 
dzz- n - x e~ z = r(-n,z), (143) 
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with 



r(-n,z) 



n—l 



3=0 



(144) 



The exponential integral function Ei(z) can be defined as 



Ei(z) 



oo /_-.\n 

7s + ln(*) + M- 

n=l 



(145) 



where = 0.5772 is Euler's constant. It reads 

T(0,z) = E 1 (z), 



(146) 



for n = 0. Close to a phase transition when the external isobaric pressure s reaches the 
hadronic bubble's internal isobaric pressure s and with the definition of Eq.(142), we find 



lim z n Y{— n, z) = 1/n, for n>0 

z— >0 

= lim (— hi(z)) , for n — 0. 

z— >0 



(147) 



Note, that the expression is logarithmically divergent for n < 0. It might be of interest for 
the reader that the integration of the delta function over the volume instead of over the 
mass will lead to the mass spectral density of the bootstrap model 



a(W, (V)) = J dvp(W, v)5 (v - (V)) . 



(148) 



On the other hand, it is wonderful to note here that it is possible to start from Hagedorn 
bootstrap mass spectral density [1, 35] to write the mass and volume density of states as 
follows 



p(m, v) = [cm~ a e bm ] 5 (m - ABv) 



(149) 



where B is the bag constant for Hagedorn bubbles. The order of the phase transition for 
the bootstrap mass spectral density with phenomenological input parameters a and b can 
be analyzed by using the isobaric partition function 



POO 

f Q (T, A; s) = C / dw-(°- 5 /2)-i e -4«(«-«o) j 

Jv 



where 



/4ET\ 3/2 



(150) 



(151) 
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The bubble's internal pressure is calculated by 

bT-1 . . 

so = B^—- (152) 

Hence, the internal MIT bag color-flavor structure parametrized with the parameter 7 is 
related to the bootstrap model by 

4 7 - 3/2 oc a - 5/2 ^7 ~ ^~^J • (153) 

The case a > 5/2 is of particular interest in the bootstrap model as it will lead (see below) 
to the phase transition. It corresponds to 7 > 3/8 in the bag model. Values of a > 7/2 
corresponds to a first order phase transition. A higher order phase transition occurs in the 
range 7/2 > a > 5/2. 

The shape and order of the phase transition for this simple model will be discussed below. 
However, we stress again that the density of states can be generalized to consider the effect 
of the bag model with a smoothed boundary where quarks and gluons are confined by the 
interacting potential in finite volume but have an extended surface. Furthermore, the above 
model can be extended, in principle, to take into account higher quantum excitations for 
Hagedorn bubbles in a straightforward fashion. 



IV. EXCLUDED VOLUME WITH SMALL AND LARGE VOLUME COMPO- 
NENTS 

The isobaric partition function for multi-particle species with a small and large volume 
components reads [36] 



[«] - ^ ex P (- [ v VdW s ]i) VsAT] m h \i 



(154) 



Z = l/ 
where 

Z = Z(T, s; mi, ^vdWi> ' ' ' > m ™' *VdW J' ( 155 ) 

and 

^exp (- KdW s L) V^K> \) = fn(T, A; {s}) + f Q (T, A; {s}), (156) 

and the set {s} = {(,h,£,q}, the parameter [s] = £h + £,q, the set { [^vdW s ] } = 
{[ v VdW s ] H ' [ v VdW s ] Q }' the Parameter [v YdW s] H = vy dWH ^ H + ^ydWtf^ and the 
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parameter [fy^w^Q = ^VcIWq^Q + v VdWQ^ H ' ^ e excmc ^ ec ^ volume effects for small and 
large components (i.e. for the known mass spectrum hadrons and Hagedorn bubbles) with 
the asymptotic approximation vq ^> v h, respectively, read 

16 , 



v V(Wh 
v VdW Q 

v VdWH 

v vaw Q 



lim 

VQ^VH 



16 3 

Q = Y 7rrQ 



[4 



1/3 + vli 3 



(V H + Vq) 
(V H + Vq) 

1/3 + t# 3 



-v H « V H 



VQ ~ VQ 



(VH + Vq) 



(157) 



The isobaric pressure for the hadronic mass spectrum reads 

= f H (T,X;{s}) 

n 

exp (- [vy dWH Z H + v VdwJo\ ) ^Stats( T ; m - A *) [ > 

n 

= [^stats( T ; m ^ A *) + < ^stats( T ; m ^^)] > 



(158) 



i=i 



where 



(159) 



Vo J mo 

oo 



are the effective excluded volume fugacity for the hadrons and their anti-particles. For the 
gas of Hagedorn bubbles, it reads 

£q = f Q (T,X,{s}) 

dvdmexp (- vy dWQ £ Q + ^y dWQ 0/ ) p(m, v)(p Q (T; m), 

dmp(m,v)(p Q (T;m) 

, 7 VdW q&3 + ^VdWg^l ) J ( u ) ' 

Vq 

POO 

/ dv e- v( ^^ H) l(v)\ . (160) 

JVo 
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= J dvex v {- [ v VdW q^Q + v VdW Q ^ H 
— dv exp ( — 

JVi 



where 



m 



x 



-B. 



(161) 



The last line in the above equation is a good approximation for sufficient large bubbles with 
Vq > m /B. The values M and V are the initial bag mass and volume, respectively. The 
fireball initial mass is taken just above the highest mass of the known hadronic mass spectrum 
particles (i.e. hadrons) listed in the particle data group book [34] M = Mq « 2.0GeV and 
the initial volume is the phenomenological volume Vq = Mq/AB where B is the MIT bag 
constant. The integration over mass density is evaluated as follows 

/•OO /"OO 

l(v) = / dmp(m,v)(p Q (T;m) = / dmp(m,v)(p Q (T;m) 

J Mr, J M n 



\ dmp(m, v) I 

Jm 

I 

J x 



mT 



M 



2tt ) 



3/2 



M 
-m/T 



dxf(x,v)e vh{x) 
where the x value is determined by 

dh(x) 
dx 



f(x,v)e vh ^K 



2ir 



—vh"(x) 



= 0. 



(162) 



(163) 



The values of f(x) and h(x) are model dependents and are determined in the Sec. V. 



A. One component excluded volume approximation 



It is worth to note that in the classical one component excluded VdW volume approxi- 
mation, the isobaric partition function reduces to 

A 1 



(164) 



E l (v?Stats(# A *) + V^StatsO 9 ; m ^ A *)) + !v dve " J ^). 
where A* = e~ VHS \i and Aj = e~ VHS X~ 1 . The isobaric pressure for Hagedorn bubbles reads 

/•OO 

f Q = / dve- vs l(v), 

I 2^ 

(165) 



/•oo 

= / dve- v[s - h{x)] f{x,' 

JVo 



-vh"(x) 
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Furthermore, in the original VdW approximation we have simply vy^-yy = vq instead of 
v VdW = ^ v °- Eq.(179) the isobaric pressure for Hagedorn bubbles in a specific model is 
given by 



POO 

f Q = C(x) / dvv-^ +l e- v[s - m] . 

The baryon density for the hadronic matter is calculated by 

hg _ T ds _n M 
Oji B n v 



where 



and 



Baryons 



+ 



T dC(x) 
C{x) dfi B 



f Q + 



T 



dh(x) 
dfi B 



C(x) / dvv-^ +2 e- v[s ~ h{x) 



V 



Baryons Mesons 



11: 



(166) 
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+ 



C(x) / (b^e- 1 ^ ) . 

Jv 



(169) 



The density for each particle reads, 



n* = n(T;mi,X*), 

/d 3 k 
T^p/statsCT;™^)- (170) 



The density of anti-particle is given by n* = n(T; m^, Aj). 



V. APPLICATIONS WITH TWO MODELS FOR THE VOLUME FLUCTUATION 

We consider the isobaric function for Hagedorn bubbles formed in the hadronic phase in 
the context of two different models. The proper choice for the bubble's volume fluctuation 
plays the crucial rule to determine the order of the phase transition for bubbles with specific 
internal color symmetries. The first model assumes the maximal volume fluctuation for 
Hagedorn bubbles in the ground state. Asymptotically, it is equivalent to the density of 
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states which is considered extensively in the literature [13-15, 23, 37]. Nevertheless, this 
choice does not lead to a real deconfinement phase transition. The second one assumes 
Gaussian-like volume fluctuation [17]. It leads to a second order phase transition to the 
deconfined quark-gluon plasma. In fact the order of the phase transition depends strongly 
on a proper choice for the volume fluctuation. The volume fluctuation is expected to be 
enhanced for low densities and high temperatures while it is supposed to be suppressed for 
high densities and low temperatures. 

A. Volume variation a la Gorenstein et. al. [13—15] 

The density of states as function of mass and volume was introduced by differentiating 
the mass spectral density (e.g. Eq.(90)) derived from the micro-canonical ensemble with 
respect to v. Nevertheless, Gorenstein et. al. [13-15] assumed that the density of states for 
a bag of unspecified volume less than v is determined by: 

a(W,(V))= f dv^-<r(W,v)= f dvp(W,v). (171) 

J(V) UV J^y) 

The volume fluctuation asymptotically behaves as the original mass spectral density 

d 

p(W, v) = ^(W, v) w P(W/v, v)a(W, v) (172) 

where 

P( X = W/V, V) ~ 1^/^-3/4^3/4 = 1 l/4 a .3/4 > ^73) 
o o 

This term is derived directly from the Eq.(172) after the terms arrangement to take the 
form of Eq.(105). The density a(W,v) is given by Eq.(90). This volume fluctuation gives 
the maximal one for < v < oo and it might be appropriate for hot and diluted matter. 
However, this case is not expected to be correct for highly compressed matter since the 
volume fluctuation is expected to be suppressed in that regime. In this case, the integral 
which appears in Eqs.(160) and (162) becomes 

T(v) = C(x)v-^ +l e vh{ *\ (174) 

where 

h{x) = \u 1/A x z / A -(x + B)/T, (175) 
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where u is given by Eq.(91) and x is the maximum point for the function h(x) and satisfies 
h{x) 



■^h(x)\ x= x = 0. It should be noted that 



x = uT 4 , 



and 



h(x) = ^Pqgp = ^ 



-uT 4 - B 
3 



(176) 



(177) 



We have also 



C(x) = ^V4^-(3 T+ l/4) ^ (x + -B)T \ 3/2 



2tt 



-h"(x)' 

3/2 



3 \ 2rr 

3 \ 2n J 

where A is given by Eq.(92). The isobaric pressure for Hagedorn bubbles becomes 

f Q (v) = £q = dve-< 4 ^)l(v), 

JVn 

-4r,+l e -v(4£ Q +€H-h(x)) ^ 



(178) 



£ Q = C{x) I dvv~ 

'Vo 

47-2 
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where 



z = V (^ Q + £ H -h(x)) 



(179) 



(180) 



Hagedorn bubble ensemble is solved self-consistently with the isobaric ensemble for the gas 
of the known hadronic mass spectrum particles 



6? — fn(T, A; £ H , £q). 



181) 



B. Volume variation d la Auberson et. al. [17] 

The basic objects are glueballs as described within the simplest version of the MIT bag 
model. In this spherical cavity approximation the energy for each glueball state is given in 
terms of the cavity volume 



rrii = 



v l/3 
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+ Bv, 



(182) 



where the bag constant energy B density simulating the confining forces is the only free 
parameter. The yi are pure numbers determined by the various modes of the eight (Abelian) 
gauge fields filling the cavity and subject to appropriate boundary conditions. The mass 
rrii and volume Vi of the standard glueball (static bags) are obtained by minimizing the 
expression with respect to the volume v. It is natural to image that the ground state 
Hagedorn bubble's volume fluctuates in the equilibrium. The mass expression is expanded 
around its minimum up to second order 

8 B 2 . o 

m ~ m H (v — Vq) . (183) 

3 mo 

The bubble's volume, v, dependence of the density of states is unknown, since the dynamics 
of the new degree of freedom v is not known. However, the mass dependence is given by the 
mass spectral density cr(W / ,f)asym and W = m — Bv. On the other hand, the function 

<p(m) ~ e- m / T , 

„ e -m 0/T . e -*§£(v-v of (lg4) 

asymptotically measures the classical volume fluctuation for the bag with mass m . The 
density of states is extracted from the statistical ensemble 



p(m ,v)ip Q (m ) ~ p(m,v)e m/T ', 



1 SB 2 ( v mp y 

cr(m — Bv, v )asyme T 3m o ^ 4B ' 



e~ m » /T . (185) 



It is worth to remind the reader the microcanonical ensemble o"(m — Bvq, v o)asym is derived 
for a bag with specific volume v and energy W = m — Bv Q . The bag's volume and mass are 
related by a strict constraint m = 4_Bf for the standard MIT bag with a sharp boundary. 
The aim of the Auberson et. al. [17] scenario is simply to soften the bag volume/mass 
constraint by allowing a small volume fluctuation m ~ 4S (vq + Sv). Therefore, we can 
write 

3 W 

p(m,v) = -P(x,v)a(x,v) aS y m , and x = — (186) 

4 v 

where 



PM = ^__ e T3m ( 4S ), 



_i — e t^+bjH 1 4s ) i 8 7) 
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where the bag energy density x = W/v is introduced. The above equation can be written 
as follows 

P(x, v) = a/ a/-7rexp [— a(v — t>o) 2 ] ~ S (v — v ) (188) 

where vq = and a = When the bag is deformed due to the highly thermal 

excitations, the smoothing function given by Eq.(112) becomes appropriate to fit the phe- 
nomenology. 

Hence with this choice we can write 



P(mMMm) = J^n- (|+|) 3/ W<.> (189) 



where 



h ^-r llv "'wr§-^ +B> ^ (190) 

with x = W/v = m/v — B. In this approximation the density of states has a narrow 
Gaussian distribution function for high temperatures and a wide one for low temperatures. 
Although this choice is suitable to describe the volume fluctuation for moderate density 
and temperature, it leads, unfortunately, to incorrect asymptotic behavior for both high 
and low temperatures. It is expected to have a very strong volume fluctuation for high 
temperature and diluted matter while the overlap effect suppresses the volume fluctuation 
for low temperature. As done in the previous section, the integration of the density of states 
over the mass is evaluated using the saddle point approximation (e.g. Eqs. (160) and (162)). 
It is given by 



X{v) = C(x)v-^ +1/2 e vh{x \ (191) 



where 



The saddle point x is the maximum point for the function h 2 (x) and satisfies -§^h(x) | x=s = 0. 
The isobaric singularity point for the fireball becomes 

£q = C(x) [°° dve-< 4 tQ+ZH-hac)) v -(4j-3/2)-^ 

JVo 

, \ 4 7 -3/2 ^oo 

= C (x) i^-j J dzz-^-^- l e-\ (193) 
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where 



zo = V (4Z Q + SH-h(x)). (194) 

C. Densities for the hadronic mass spectrum particle gas and the Hagedorn bubble 

gas 

The pressure for a gas comprising the mass spectrum of all known hadrons and fireballs 
reads 

P = T(Z H + Z Q ), (195) 

Here the fireballs are simply the Hagedorn states. The total isobaric pressure is calcu- 
lated from the extreme right singularity as given Eqs.(132) and (133) and demonstrated by 
Eq.(154) for small and large volume components. The isobaric pressures for the gas of the 
hadronic mass spectrum particles and the gas of the fireballs are given, respectively, by 

i 

£ Q = [°° dve-< 4 ^ H )l(v), (196) 

Jv 

where A* = e ~ VH \^ H+ ^Q) \ and \ = exp (Bifi B /T + Sifis/T + •••). The baryonic density 
for the hadronic gas can be calculated by differentiating the total pressure with respect to 
the baryonic chemical potential 

„ _ Spectrum .Fireballs 
n B — n B -+- n B 

= T-^-tn + T-^-tQ. (197) 

VI. THE ORDER OF THE PHASE TRANSITION TO QUARK-GLUON 
DROPLETS OR PLASMA STATE 

The order of the phase transition from the hadronic phase to quark-gluon droplets or 
plasma state is determined by scrutinizing the properties of the isobaric pressure for Hage- 
dorn bubbles when integrating over the volume. The volume fluctuation in the model is cal- 
culated after rescaling the mass and volume (W, v) variables in the (grand-) micro-canonical 
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density to mass density and volume (x = W/v,v) variables. This mass density/ volume 
scaling is justified by the assumption that v oc W and v — ^| for the MIT bag mode. 

Although Hagedorn bubble's internal color-flavor symmetry is important to determine 
the order and shape of the phase transition, it is not the only the criteria. The bubble 
volume fluctuation after mass/volume scaling is, indeed, crucial to fix the order of the phase 
transition, in particular for hot and diluted hadronic matter. It is reasonable to expect that 
the volume fluctuation varies differently for compressed matter where the bubbles start to 
overlap each other, and for diluted and hot matter where the bag's surface is expected to 
dissociate spontaneously near the critical temperature. The hadronic phase comprises of all 
the known hadronic states (mass spectra of baryons, mesons and their resonances) as well as 
the highly excited hadronic fireballs, i.e. Hagedorn bubbles with higher internal color-flavor 
symmetry. Each hadronic fireball is approximated as an ideal gas of quarks and gluons 
moving freely inside the bag within a specific color-flavor quantum state. 

However, our numerical calculations show that Hagedorn bubbles appear only in highly 
dense matter for large baryo-chemical potential. In contrast they are unlikely to appear in 
diluted matter for low baryo-chemical potential even when the system is approaching the 
critical temperature. This observation apparently contradicts to the common thought that 
Hagedorn bubbles always show up below T c for dilute and hot hadronic matter. Usually 
Hagedorn states are supposed to develop below T c . However this apparent contradiction is 
resolved easily by noting that the bubble's internal pressure is quartic temperature dependent 



This pressure grows up and exceeds the external hadronic pressure quickly with respect to 
the temperature. Below the critical temperature, these Hagedorn bubbles are suppressed 
by the external hadronic pressure of the mass spectrum gas. Furthermore, in the numerical 
calculations we have taken a relatively small Van der Waals repulsive volume for the mass 



suppresses these Hagedorn bubbles with relatively large excluded volumes strongly. The nu- 
merical calculations also show that smaller values for the initial Hagedorn bubble's excluded 
volume enhance the appearance of Hagedorn states significantly at temperatures below T c 
for diluted hadronic matter. The small excluded volume for the mass spectrum particles 




(198) 




■ = 4 x 0.0654fm 3 ) while a relatively large one for Hage- 
4 x 4.190fm 3 ). This small-large excluded volume scenario 
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and the large excluded volume for Hagedorn states effect concerns us in the present work in 
order to demonstrate clearly the different scenarios for the phase transition diagrams. These 
relatively large Hagedorn bubbles can in principle appear in the system due to the high ther- 
mal excitations of the vacuum. Whenever they appear in dilute hot matter, however, they 
are mechanically unstable. The bubbles with an internal pressure less than the external one 
collapse and disappear rapidly. However, when the internal pressure of Hagedorn bubbles 
reach the external one from below near the critical temperature, they will eventually expand 
and explode because the bubble's overlap affect is negligible in this regime. The resulting 
exploding bubbles expand rapidly forming quark-gluon droplets by merging with each other 
eventually filling the whole space. This big quark-gluon droplet loses its internal structure 
as it expands and undergoes a true deconfinement phase transition. On the other hand, 
Hagedorn bubbles can exist in the compressed and cold hadronic phase. Their appearance 
at large baryonic chemical potential is essential to soften the equation of state being a mea- 
sure of the highly excited mass spectrum of compressed hadronic matter. These bubbles 
expand slowly in dense matter because of the overlap with other bubbles. As the system is 
compressed further at low temperature, many bubbles are likely to merge into each other 
to form denser Hagedorn bubbles with more complicated internal color-flavor structure. As 
this effect increases with baryon density, the bubble volume fluctuation will decrease corre- 
spondingly. Therefore, it is expected that the hadronic system undergoes a phase transition 
to foam of highly dense Hagedorn bubbles at large baryonic chemical potential and low 
temperatures. Furthermore, when the foam is heated due to e.g. compression, the internal 
surfaces dissolve the bubbles merge to form bigger droplets. At some point the quark-gluon 
droplet surface collapses and the system undergoes a phase transition to fully deconfined 
quark-gluon plasma. 

The integration over volume for Hagedorn bubble's isobaric pressure reads 



where Zq = Vo(4£q + £# — s )- The exponent a is some phenomenological parameter depend- 
ing on both the bubble's internal color symmetry and the volume fluctuation of the bubble 




(199) 
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itself. The asymptotic behavior near the phase transition point acts as 

lim S(— a,zo) ~ lim z^^Ti— a + 1, zq) — > finite for ct > 1 

2„^0+ 2 ->0+ 

~ lim ln(zo) — > diverges for a = 1, 

20— >0+ 

— > diverges for a < 1, (200) 

where z = V^(4£q + — s ). The factor 4 which appears in front of £q comes from the hard 
core Van der Waals repulsion for large and small components for Hagedorn bubbles and the 
mass spectrum particles, respectively, as discussed in Sec. IV. Therefore, the convergence or 
divergence of the isobaric pressure near the point of the phase transition can be summarized 
as follows 

lim £q — > diverges (a — 1 < 0), 

2 ^0 + 

lim £ Q -> finite {a - 1 > 0). (201) 

Note, that the phase transition does not exist for Hagedorn bubbles characterized by volume 
structure a < 1. Hagedorn bubble's external pressure diverges and subsequently the internal 
pressure is always less than the external one and consequently these bubbles collapse and 
are strongly suppressed in the hadronic phase. 

The phase transition occurs only for bubbles with a volume parameter a > 1. Nonethe- 
less, one can expect a rapid and smooth phase transition for bubbles with internal structure 
a = 1 + . In this case the bubble external pressure is large but finite near the phase transition. 
Hence, it is possible that expanding Hagedorn bubbles emerge and form big quark-gluon 
droplet which occupies most of the available space. The instability of Hadronic bubbles 
increases when the hadronic external pressure increases to large values just below the phase 
transition. The phase transition will take place when Hagedorn bubble's internal pressure 
becomes equal to the external pressure of the hadronic gas 

^ q + £h = s . (202) 

Despite the quark-gluon droplet forming rapidly and expanding quickly, not the whole 
hadronic gas undergoes a phase transition. In these circumstances, it would be difficult 
to distinguish between the hadronic phase and the quark-gluon plasma and the system un- 
dergoes a smooth but rapid cross-over phase transition. The order of the phase transition is 
determined by the discontinuity of the n-th derivative of the isobaric pressure at the surface 
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of Hagedorn bubble just below the phase transition line. If the first derivative is discontin- 
uous, then the system undergoes a first order phase transition. When the first derivative is 
continuous, then the second derivative should be inspected. If the first derivative is contin- 
uous while the second derivative is discontinuous, then the system undergoes a second order 
phase transition. Furthermore, when the first and second derivatives are both continuous 
while the third derivative is discontinuous then a third order phase transition takes place 
and so on. Consequently, the n th order phase transition is determined by the discontinuity 
of the n th derivative of the isobaric pressure at the surface of Hagedorn bubble. The first 
derivative 

(4£q + in ~ so)' oc ^/S(-a + 1, z ) oc l/z^ 2 T(-a + 2, z ), 

oc *o~ (a ~ 1)+ \ *o = V (4£q + £ h - s ) , (203) 

at the surface of Hagedorn bubble is continuous only with the volume parameter — (a — 
1) + 1 > 0. The prime notation in X' indicate the partial derivative with respect to the 
thermodynamical ensemble such as the temperature and chemical potential (e.g. X' 
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The continuity of the second and third derivatives 

(4£q + £h~ s )" oc E(-a + 2, z )/E 3 {-a + 1, z ), 

oc ^" 2(a - 1)+1 ,z = Vo (4£q +Zh-s ), (204) 

and 

(4£q + £h - so)'" oc [E(-a + 3, z )E(-a + 1, z ) - 3E 2 (-a + 2, z )] /S 5 (-a + 1, z ), 

oc z~ 2{a - 1)+ \ z = V (4^q + £ h - so) , (205) 

appears only for — 2(a — 1) + 1 > and — 3(« — 1) + 1 > 0, respectively. The discontinuity 
of the n th derivatives 

(4£q +£ h - s ) (n) oc lim \n(z ) -> diverges, z = V (4^ Q + ^ - s ) , (206) 

are determined by — n(a — 1) + 1 < 0. The n th derivative is continuous 

(4£q + £ h - so) {n) oc lim z " n(a - 1)+1 ~ 0, (207) 

zo-*0 
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for — n(a — 1) + 1 > 0. Hence the condition for the n-th order phase transition reads 

a > 1 + l/n. (208) 

The 1 st , 2 nd , n th order phase transitions are given by«>2, 2>a>l + | and [1 + ^-j-] > 
a > [l + i] , respectively. However, a cross-over phase transition takes place only when 
all derivatives of n th order are becoming equal. This corresponds to the volume parameter 
a = 1 + . The external pressure of Hagedorn bubble diverges smoothly as a approaches 1 + 
and the system subsequently undergoes a phase transition. 

The order of the phase transition depends on both the bubble volume fluctuation and its 
internal symmetry. We demonstrate this dependence by considering the two different models 
for the volume fluctuation as presented in section III. The volume fluctuation structure as 
given by Gorenstein et. al. [13-15] (see Eq.(172)) with specific internal color symmetry is 
handled by the redefinition 

a = 4 7 -l. (209) 

In this approach, the ground state Hagedorn bubble expands freely to infinity, i.e. the 
maximal volume fluctuation for the bubble in the ground state. We have find then that the 
phase transition is not possible for colored bag of gluons and flavorless quarks with structure 
7 = \ as a — 1. The first and second order phase transitions are possible for bags with 
color structures characterized by 7 > | and 7 > |, respectively. Furthermore, the n th order 
phase transition takes place at 7 > The colorless bag of gluons and flavorless quarks 

has 7 = | undergoes a first order phase transition. 

On the other hand, the Gaussian volume fluctuation given by Eq.(186) and suggested by 
Auberson et. al. [17] corresponds to 

a = 4 7 - 1/2. (210) 

for Hagedorn bubbles with specific internal color-flavor symmetry. Hence, the bubbles with 
color-flavor structure 7 = 1/2 undergo a third order phase transition to the real deconfined 
quark-gluon plasma as (a = 3/2). The second order phase transition occurs for bubbles with 
§ > 7 > I while those with 7 > | have a first order phase transition. Hence, the Hagedorn 
bubbles in the color singlet state and internal color-flavor structure 7 = 3/2 (a = 11/2) 
undergoes a first order phase transition. 
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Generally speaking, the volume fluctuation basically depends on the quantum wave- 
function of the quark-gluon cluster and varies with respect to temperature and chemical 
potential. However, for the bag model with a deformed boundary, the bag volume is ex- 
pected to fluctuate around the mean volume (v — v o) uniformly which can be approximated 
by a Gaussian for some range of chemical potentials. Nonetheless, this approximation is not 
necessarily correct for the entire \i — T phase diagram. The Gaussian approximation fails 
to measure the real volume fluctuation for Hagedorn bubbles embedded in the compressed 
matter in particular when these bubbles overlap with each other. In this regime, the bub- 
bles squeeze into each other and there is a little room for further expansion. Hence, the 
volume fluctuation is likely damped for highly compressed and cold matter. Moreover, the 
Gaussian approximation also fails for hot and dilute baryonic matter regime as the hadronic 
bag surface dissociates spontaneously near the critical temperature and Hagedorn bubbles 
expand rapidly and resulting in a first or even higher order phase transition to quark-gluon 
droplets or real deconfined quark-gluon plasma. 



VII. RESULTS AND DISCUSSIONS 



In the following, the hadronic matter is treated as Van der Waals gas consisting of all 
the known particles as given in the particle data group book [34]. Hereinafter, we call these 
particles the mass spectrum particles. The highly excited hadronic states are also taken 
into account and they are taken as Van der Waals gas of Hagedorn bubbles. The density of 
states for Hagedorn bubbles is derived from Laplace inverse of the mixed canonical ensemble 
of blob of quarks and gluons with specific internal color-flavor structure. These states are 
similar to Hagedorn states in the bootstrap model [1, 35] where the density of states for the 
hadronic matter is modified significantly for the highly compressed or heated matter because 
the appearance of exotic hadronic (i.e. Hagedorn) bubbles with large masses m# > 2.0 GeV 
in the hadronic phase. The rich mass spectrum of Hagedorn states is essential for softening 
the equation of state in particular for a large baryonic chemical potential. In order to 
demonstrate the importance of Hagedorn bubbles in the hadronic phase, we have studied 
the quark-gluon plasma phase transition with and without Hagedorn bubbles. The effect of 
the initial volume fluctuation of Hagedorn bubbles in the phase transition is also considered. 
The density of states for Hagedorn bubbles is considered for the quark and gluon bags with 
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various internal color-flavor symmetries range from colored bubbles (or colorless bubbles with 
the minimum color-flavor correlation) to singlet ones with strong color-flavor correlations 
specified by the parameter 7 co i or < 7 < 7 s i n glet ^ ^strong an< ^ * n the context of two 
approaches for the bubble's volume fluctuation. The first approach is based on Gorenstein 
et. al. [13-15] ansatz for the maximum volume fluctuation while second one is based on 
Auberson et. al. [17] ansatz for Gaussian like volume fluctuation. The bag constant for the 
hadronic pressure is taken B 1 ^ = 210MeV and the initial bag volume starts at |7r (l.Ofm) 3 . 
In Figs. 1, 2 and 3, the one exclude volume component approximation is considered in the 
numerical calculations while the excluded volume with small and large volume components 
approximation is considered in Figs. 4, 5, 7 and 8. 

We display the baryonic density and temperature (p# — T) phase transition diagram in 
Fig. 1. The thin lines depict the phase transition diagram for hadronic matter consisting of 
Van der Waals gas of the known particles without Hagedorn bubbles with various excluded 
volumes. In this approximation, the phase transition to quark-gluon plasma is calculated 
using the standard Gibbs construction scheme where the pressures and the chemical poten- 
tials are equal in both phases. The quark-gluon plasma is treated as an ideal gas of quarks 
and gluons. The excluded volume for the mass spectrum particles is taken ^vdWn = ^ Vr 
where Vh = and the phase transition diagrams are displayed for excluded volume radii 

rjj=0.25, 0.35, 0.40 and 0.45 fm. The number of particles per volume is reduced as the size 
of the excluded volume increases. The proper Van der Waals excluded volume is chosen to 
fit the phenomenology and it is taken to be rather small ^vdWn = ^^ r H ~ 0.40-0.65 fm 
[r H ~ 0.25-0.45 fm). The classical Van der Waals excluded volume corresponds the phe- 
nomenological nucleon radius 0.40-0.65 fm. On the other hand, the thick lines display 
the baryonic density for the hadronic matter consisting all the known hadronic spectrum 
particles in the particle data group book [34] as well as the highly excited Hagedorn bub- 
bles. The excluded volume for the known hadronic particles is taken to be r H = 0.25fm, i.e. 
r VdW = 0-40fm. The density of states for Hagedorn bubbles is taken as given in Eqs(90), 
(172) and (173). The initial volume fluctuation, Vo, for Hagedorn bubbles is taken rpg=1.0 
and the bubble's Van der Waals excluded volume VyjWg — 4Vo- These Hagedorn bubbles 
are also supposed to be the fireballs (FB) could appear in the hadronic phase. The bag 
constant for the Hagedorn bubble is taken B 1 / 4 = 210MeV. The internal symmetry for 
the color-flavor correlation in Hagedorn bag is chosen to be one for the massless flavorless 
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and color singlet state, 7=1.5 with a smoothed volume fluctuation as given in Eq.(173). 
When Hagedorn bubbles are included in the calculation, the phase transition diagram is 
found self-consistently by searching the second singularity in the isobaric pressure. The 
first singularity is the extreme right singularity point in the limit of infinite external vol- 
ume given by Eqs.(132) and (133), while the second singularity is for the isobaric pressure 
of Hagedorn bubbles. This second singularity is analogous to the condition for the phase 
transition determined by Hagedorn bubble's mechanical instability arises when the bubble's 
internal pressure reaches the external one from below. Due to this pressure instability, the 
bubbles start to expand rapidly and fill the entire space, consequently, quarks and gluons 
move freely and the phase transition to the quark-gluon plasma phase is reached. Fig. 1 
shows that although the Van der Waals gas of the known hadronic mass spectrum particles 
without Hagedorn bubbles gives a reasonable phase transition diagram for low and inter- 
mediate densities, it fails to predict the phase transition for large baryonic density and low 
temperature. The finite size effect for the hadronic phase and quark-gluon plasma is also 
found important here [38-40]. Including Van der Waals gas for color singlet bubbles 7 = § 
and bubbles with more complicated internal color-flavor structure 7 = 3.0 shifts the phase 
transition line to higher baryonic densities. The phase transition diagram exhibits a long 
"tail" at low temperatures and high densities. Furthermore, it is found that increasing the 
bubble's initial volume suppresses Hagedorn bubbles population. This supports the intuitive 
idea that Hagedorn bubble spectrum is the continuation for the excited mass spectrum for 
the known hadronic particles found in the particle data group book [34] where Hagedorn 
bubble mass spectrum starts at the end of the hadronic mass spectrum particles. 

The above thin lines show the baryonic densities for the quark-gluon plasma phase p% GP 
above the phase transition point for hadronic matter consisting both the spectrum of known 
hadrons and Hagedorn bubbles with various internal color-flavor structures. It is shown 
that Pg GP above the phase transition point exceeds significantly that for the hadronic gas 
Pg G consisting the known hadronic particles and Hagedorn bubbles with internal structure 
7 = 1.5 below the phase transition point. This measurable discrepancy in the baryonic 
density between two phases indicates a first order phase transition and a discontinuity in 
the baryonic density at the point of phase transition. However, the case is rather different for 
hadronic matter consisting Hagedorn bubbles with internal color-flavor structure 7 = 0.51 
0.5833 and 0.625. It is found that the split in baryonic density between two phases is 
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tiny for 7 = 0.51 and small for 0.625. This could indicate a continuity for the baryonic 
density at the point of the phase transition and a discontinuity for higher derivatives of the 
thermodynamical grand potential density. Indeed the baryonic density split above and below 
the phase transition line sheds some information about the order of the phase transition 
and the discontinuity of the isobaric pressure and its first n-th derivatives. Indeed, it is 
supposed that the higher order phase transition reduces the split size at the point of the 
phase transition significantly where the isobaric pressure is continuous while its n th derivative 
is supposed to be discontinuous. Furthermore, the numerical calculations show that this 
density split decreases significantly as the internal structure decreases from 7 = 1.5 to 0.51. 
This evidently indicates that the hadronic system consisting Hagedorn bubbles with internal 
structure 7 = 1.5 undergoes a first order phase transition while the system consisting bubbles 
with internal structure 7 = 0.51 undergoes a higher order phase transition and the n th higher 
derivatives of the thermodynamical grand potential density is continuous. 

The phase transition diagram in the baryonic chemical potential and temperature (T—/ib) 
plane is displayed in Fig. 2. The phase transition line for the hadronic gas of only the 
known hadronic mass spectrum particles without any Hagedorn state ends at the chemical 
potential pb ~ 1300 MeV for low temperatures. However, the size effect for the excluded 
volume in T — ps plane is not apparent as for the (p# — T) phase transition diagram shown 
in Fig. 1. When Hagedorn bubbles are included, the phase transition diagram is shifted to 
larger baryonic chemical potentials for low temperatures. The phase transition diagram is 
modified to have a long tail at large baryonic chemical potentials and low temperatures. 
When Hagedorn bubbles appear in highly dense matter and low temperatures, the system 
prefers to stay in the hadronic phase which is dominated by highly compressed Hagedorn 
bubbles. Furthermore, the system experiences a smooth phase transition from hadronic 
gas dominated by the hadronic mass spectrum particles to another one dominated by foam 
of dense Hagedorn bubbles. The baryonic densities for both the mass spectrum particles 
pHG.spectrum jjagedorn bubbles p^ G ' FB in the hadronic phase below the phase transition 
line to QGP are displayed in Fig. 3. We have also displayed the baryonic density p% GP for 
QGP above the line of the phase transition. It is shown that the hadronic gas of the hadronic 
mass spectrum particles dominates the hadronic phase for low and intermediate baryonic 
chemical potentials with temperatures higher than 60 MeV while Hagedorn bubbles become 
the dominant one in the hadronic phase for the chemical potential ps exceeds 1200 MeV 
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and the temperature falls below 60 MeV. 

As outlined in the previous sections, the density of states for Hagedorn bubbles is given 
by the microcanonical ensemble for blobs of quarks and gluons with specific color-flavor 
symmetry. The micro-canonical ensemble is derived from the inverse Laplace transform of 
the mixed grand canonical ensemble of quarks and gluons confined in a cavity. The micro- 
canonical ensemble measures the mass spectral density Eqs.(80) and (90) for Hagedorn 
bubble with specific volume for a cavity with sharp surface. The volume fluctuation can be 
calculated by finding a solution for a bag model with a deformed boundary and this solution 
is, in general, not known. The bag model with a smoothed boundary can be mimicked by 
using several assumptions to fit the phenomenology. For the nuclear shell model, the density 
of states is modified by using the ansatz of Strutinsky [27-31] in order to fit the nuclear data 
and is based on smearing the delta function 5(E — E ) respecting energy conservation. 

Here, we are going to study different scenarios for incorporating volume fluctuation and 
the phase transition impact on extending our previous discussions to more realistic systems. 
We emphasize in the following that the bubble's color-flavor internal symmetry and its 
volume fluctuation plays a vital role in determining the order of the phase transition at low 
baryonic chemical potential and high temperature and the shape of the phase transition line 
at large chemical potential as well. In order to emphasize the volume fluctuation's role in 
the phase transition we studied two models with different ansatz for the volume fluctuation. 

In the model of Gorenstein et. al. [13-15] hereafter denoted as model (I), the volume 
fluctuation is measured by differentiating the microcanonical ensemble with respect to the 
bag's volume (see Eq.(171)). This approximation leads to a smoothed volume fluctuation 
as given in Eq.(173) which is independent on the bubble's volume and depends only on 
the bubble energy density. It is asymptotically equivalent to replacing the volume delta 
function S(v — v ) given in Eq.(137) in the standard MIT bag by the power law x 3//4 where 
x is the energy density. This variation is volume independent and subsequently it behaves 
the maximum volume fluctuation. We display the phase transition diagram in Fig. 4 with 
various values for Hagedorn bubble's color-flavor correlation parameter, 7= 0.51, 0.5833, 
0.625, 0.75, 1.0, 1.5, 2.0 and 3.0. The excluded volume for the mass spectrum particles is 
fixed to rjj = 0.30fm and Hagedorn bubbles' initial volume fluctuation starts from rpg = 
l.Ofm. The 7 = 0.5 case corresponds to colored bubbles while 7 = 1.5 corresponds to color 
singlet bubbles. For values of 7 > 1.5 Hagedorn bubbles have more complicated color-flavor 
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structures but still remain in a color singlet state, while for 7 < 1.5 the bubbles have color- 
flavor structure with color nonsinglet components. The cases with 7 = 0.51, 0.5833, 0.625, 
0.75 or 1.0 do not alter the shape of the phase transition line for large baryonic chemical 
potential. The phase transition line ends at a baryonic chemical potential of //# ~ 1350 MeV 
for zero temperature. The order of the phase transition changes from 1st order to 2 nd , 
3 rd , A th and n th order for 7 = 0.75, 0.625, 0.5833 and 0.51, respectively. For 7 = 1.5, 
2.0 and 3.0, the phase transition diagram changes significantly and the phase transition 
line is shifted to larger baryonic chemical potential for low temperature. The shift of the 
phase transition line increases drastically with 7. However, in the one component excluded 
volume approximation (see Fig. 2) this shift in the phase transition diagram to larger baryonic 
chemical potential at low temperature is less pronounced for 7 = 1.5 than that for the small 
and large excluded volume components approximation displayed in Fig. 4. The reason is 
that the bubbles in the one excluded volume component approximation are effectively more 
suppressed than that for the small and large excluded volume components approximation. 
The system prefers to remain in the hadronic phase dominated by Hagedorn bubbles for 
large baryonic chemical potential rather than undergoes a phase transition to real deconfined 
quark-gluon plasma. When the medium becomes sufficiently hot, these Hagedorn bubbles 
undergo a phase transition to the quark-gluon plasma. 

We display in Fig. 5 the pressures for the gas of all known nonstrange particles and the 
gas of Hagedorn bubbles versus the baryonic chemical potential \ib at temperatures just 
below the phase transition diagram. The sum of both pressures gives the total pressure for 
the hadronic phase. Most of the thermodynamical quantities such as particle multiplicities 
are derived from the isobaric pressure. The thermodynamical grand potential density is 
also given in terms of the isobaric pressure. In the dilute hadronic matter, the pressure for 
Hagedorn bubble gas is suppressed and Hagedorn states unlikely appear in the hadronic 
phase. The gas of the hadronic mass spectrum particles is the dominant in the hadronic 
system for baryonic chemical potential up to \xb ~ 1200 MeV. In this regime the colored 
bubbles with 7 = 0.5 are strongly suppressed. If these bubbles appear because of the high 
thermal excitations, their gas pressure diverges and whenever they appear they explode 
and occupy the entire space and mix with the hadronic mass spectrum particles instantly. 
This mechanism leads to the so called crossover phase transition. The bubbles with color 
nonsinglet components and exponents 7 < f may appear slightly just below the phase 
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transition line and the system undergoes a higher order phase transition to quark-gluon 
plasma. These bubbles are also suppressed in the hadronic phase, their appearance just 
below the phase transition diagram allows a higher order phase transition. On the other 
hand, the bubbles with the net color singlet state do not appear for small and intermediate 
baryonic chemical potentials < 1000 MeV. When the baryonic chemical potential exceeds 
fiB > 1000 MeV, these Hagedorn bubbles appear and dominate the hadronic phase. For the 
small and intermediate baryonic chemical potential, Hagedorn states appear just below the 
phase transition line. Fig. 5 shows also that with increasing 7, the pressure for Hagedorn 
bubbles increases and the system prefers to remain in the hadronic phase for dense and cold 
nuclear matter. However, at high temperatures, the system undergoes a phase transition to 
quark-gluon plasma. 

In order to summarize the scenario for the order and shape of the phase transition in 
model (I), the corresponding phase transition diagram is sketched in Fig. 6. Let us discuss 
the influence of the parameter a = 47 — 1 on changes to the phase transition diagram 
by varying it with the baryochemical potential. There is no phase transition for hadronic 
matter with a = 1 (e.g. bubbles structure 7 = 0.5) and the system undergoes a crossover 
phase transition to the quark-gluon plasma for dilute and hot matter. For slightly larger 
baryonic chemical potential, the system undergoes n-th order phase transition to quark-gluon 
plasma for a < 2 (7 < |) as Hagedorn bubbles appear near the phase transition due to the 
thermal excitations. The first order phase transition takes place only for bubbles with a > 2. 
Hadronic matter consisting bubbles with a = 2 undergoes a second order phase transition 
while the hadronic matter consisting ones with a = 3/2 (7 = 0.625) undergoes a third order 
phase transition. Therefore, the hadronic matter consisting color singlet bubbles, which 
has internal structure 7 = 3/2 for gluons and massless flavorless quarks undergoes a first 
order phase transition. For example, the bubbles with specific baryonic and mesonic internal 
structures have 7 = 3/2 and 7 = 13/8 [18], respectively, and bubbles with 6 massless flavors 
and 3 colors with exact SU Gjt —> SU(3) x SU(3) have 7 = 6 [18]. Therefore, it is expected 
that the strong color-flavor correlation will lead 7 ^> 1 for three flavors with symmetry 
SU(3) C x SU(3) fL x SU(3) fR and for two flavors with symmetry SU(3) C x SU{2) fL x SU(2) fR . 

The basic question is we follow now at which extent do the color-flavor symmetry and the 
modification of the volume fluctuation affects the order and shape of the phase transition 
line. It is known from the nuclear shell theory [27-31] that the deformation of the cavity 
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boundary smears the density of states by smoothing the boundary surface. Auberson et. 
al. [17] have argued that within the MIT bag model the mass fluctuates around the initial 
bag volume and it is approximated by Eq.(183). This procedure leads to the Gaussian-like 
volume fluctuation (see Eq.(187)). The resulting density of states for Hagedorn bubbles is 
given in Eq.(186). The quark and gluon bubble's Gaussian volume fluctuation approach will 
be denoted as model (II). The phase transition diagram in the (T — /ib) plane is displayed in 
Fig. 7 for various internal color-flavor structures 7 = 0.5, 1.0, 1.5, 2.0 and 3.0. The excluded 
volume for the hadronic mass spectrum particles is taken to be = 0.30fm. The bubble's 
initial volume fluctuation starts from rpg = l.Ofm. 

Fig. 7 exhibits a first order phase transition for the colorless quark and gluon bubbles with 
7 = | and also for bubbles with color nonsinglet components of | > 7 > |, a second order 
phase transition for bubbles | > 7 > \. A third order phase transition is found for colored 
bubbles with 7 = |. This means that the third order phase transition is the maximum 
order. It takes place near the critical temperature T 155 MeV where the hadronic mass 
spectrum gas undergoes a third order phase transition to real deconfmed quark-gluon plasma. 
Hagedorn bubbles with a strong color-flavor correlation 7 > | but with net color singlet 
states trigger always a first order phase transition. Furthermore, the phase transition line 
for bubbles with color nonsinglet components and 7 < | are almost the same, although they 
undergo phase transitions of different order. The phase transition line drops at a chemical 
potential of /is ~ 1350 MeV for low temperatures, T < 100 MeV. The interesting result to 
note in this scenario is, that the bubbles with specific internal color-flavor correlations and 
specific color components undergo a second order phase transition and that only colored 
bubbles undergo a third order phase transition. We would like to remind the reader that a 
specific flavor correlation has been introduced via the center of mass projection in order to 
reproduce the bootstrap density of states [24]. The Lorenz structure of the center of mass 
projection is not invariant in the hot bath and can be broken partially or completely. This 
mechanism might modify the density of states. In the present model, we have assumed that 
this mechanism, in somehow, is absorbed by the phenomenological color-flavor parameter 7. 
The story is rather different for color singlet quark and gluon bags (i.e. Hagedorn bubbles) 
but with complicated color-flavor internal symmetries. Hagedorn bubbles with 7 = 1.5, 2.0 
and 3.0 modify the phase transition diagram significantly and the system prefers to stay in 
the hadronic phase for large chemical potential and sufficiently low temperature. When the 
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medium's temperature exceeds T = 100 MeV, the bubble color-flavor symmetry becomes 
less important to modify the shape of the phase transition line although it plays the vital role 
in determining the order of the phase transition. Hence, the phase transition diagram appear 
to be almost the same for Hagedorn bubbles with various color-flavor internal symmetries. 
This means that regardless of the bubbles internal color-flavor complications, quark-gluon 
droplets or plasma is formed. 

We display in Fig. 8 the pressures for the gas of known spectrum particles and the gas 
of Hagedorn bubbles versus the baryonic chemical potential //# at temperatures just below 
the phase transition line. Hagedorn bubbles are strongly suppressed in diluted hadronic 
matter and they don't appear at small baryonic chemical potential. They likely appear 
just below the phase transition line because of high thermal excitations of the medium. 
Furthermore, when the baryonic chemical potential increases and exceeds \ib ~ 1000 MeV, 
Hagedorn bubbles start to appear in the hadronic phase and their gas pressure increases 
significantly as the baryonic chemical potential increases. They becomes more dominant 
than the hadronic mass spectrum particle gas for highly compressed hadronic matter. 

The sketch of the phase transition diagram for model (II) is depicted in Fig. 9. The 
parameter a = (47 — |) is assumed to be modified in the medium and increases as the 
hadronic matter is compressed. In this scenario, the hadronic phase undergoes a third order 
phase transition to colored bubbles with a — 1.5 (7 = 0.5) and forms real deconfined quark- 
gluon plasma in diluted and hot matter. The bubbles with color nonsinglet components 
with structure 2 > a > | (i.e. (| > 7 > |)) have a second order phase transition while 
for y > a > 2 ((§ > 7 > §)) there appears a first order phase transition. The hadronic 
phase consisting Hagedorn states given by color singlet bubbles (7 > 3/2) undergoes a first 
order phase transition to quark-gluon droplets or plasma. Increasing the value of 7, modifies 
the phase transition diagram shape substantially in particular when the baryonic chemical 
potential becomes sufficiently large. 

The comparison between phase transition scenarios for models (I) and (II) shows that 
the modification of Hagedorn bubble's internal structure due to the interquark potential is 
important in order to draw and determine the order of the phase transition diagram. The 
interquarks potential modification could slightly deform the bag boundary and consequently 
the order of the phase transition even if the internal color-flavor symmetry remains intact. 

Naively, Hagedorn bubble's internal structure a (or 7) can also be modified in the 
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medium. The rough sketch for the dependence of the phenomenological Hagedorn bub- 
ble's internal structure 7 on the volume, temperature and chemical potential is depicted in 
Fig. 10. Primarily insight on melting the frozen internal color degrees of freedom in the 
medium have been given by Elze and Greiner [20]. 

It is possible to go beyond the approximations of the previous models by considering 
higher order volume fluctuations or even finding a realistic hadronic bubble's wavefunction. 
In this context, it is reasonable to assume that Hagedorn bubbles are thermally excited to 
higher quantum states whenever they appear. These highly excited Hagedorn bubbles can 
then even evaporate or emit smaller hadronic bubbles. In this case, the volume fluctuation 
becomes stronger and the bubble's internal color-flavor structure is modified by smoothing 
the cavity boundary and consequently the density of states for Hagedorn states. When the 
phase transition takes place, the quark-gluon droplets continue to expand in the hot and 
diluted matter and then they eventually overlap with each other to form real deconfined 
quark-gluon plasma. On the other hand, in dense and cold hadronic matter the situation 
is rather different, however. When the baryonic chemical potential increases, the hadronic 
bubbles agglomerate and merge to form bubbles which are likely to have higher color-flavor 
symmetry and larger baryonic number. Hagedorn bubbles with simple symmetries merge 
and form dense bubbles with higher color-flavor symmetries and soften the equation of 
state. They shift the phase transition to higher density and temperature. At warm and large 
baryonic density, the system is dominated by Hagedorn bubbles. When these dense Hagedorn 
bubbles expand and overlap with each other, the system will undergo a phase transition to 
quark-gluon droplets. These droplets continue to expand and lose their internal color-flavor 
symmetry and then eventually they merge all together to form quark-gluon plasma. 

The baryonic number for each Hagedorn bubble increases when the system is compressed 
and cooled and more dense Hagedorn bubbles with larger baryonic numbers appear in the 
system. At highly compressed matter, most of the baryonic density tends to concentrate 
in Hagedorn bubbles. However, the bubble's size fluctuates smoothly and even shrinks in 
order to reduce the overlap effect with other hadronic bubbles. Hagedorn Bubbles with low 
mass density merge with other bubbles to form denser ones. When the maximum density 
is reached, the system undergoes a phase transition to Hagedorn bubbles foam. For low 
temperature and very dense matter, the bubble's size can not shrink anymore due to the high 
constituent quarks and gluons pressure and the bubbles high surface tension in particular 
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for cold matter. Hence the equation of state can be softened only by merging Hagedorn 
bubbles to form denser bubbles. Furthermore, the thermal excitation dissolves the surface 
between the bubbles spontaneously and subsequently the Hagedorn bubbles foam collapses 
to form quark-gluon plasma. 

VIII. SUMMARY AND CONCLUSIONS 

We have studied the order, shape and critical point for the phase transition diagram for 
hadronic matter consisting of all the known nonstrange hadronic particles and the highly ex- 
cited hadronic bubbles correspond Hagedorn states those existed in the extreme conditions. 
The basic assumption is that Hagedorn bubble's internal structure exponent 7 depends on 
the medium and may modify itself self-consistently with respect to /is in particular for tem- 
perature just below T c . The exponent 7 is related to Hagedorn bubble's internal color-flavor 
potential interaction. We have shown the order of the phase transition depends basically on 
the phenomenological exponent 7. We have demonstrated that both the volume fluctuation 
and the internal color-flavor structure for Hagedorn states play a crucial role in determining 
the order of the phase transition to quark-gluon droplets or plasma. It is found that Hage- 
dorn bubbles' quantum excitations modify the volume fluctuation by smoothing the quarks 
and gluons cavity boundary and subsequently change the order of the phase transition. The 
excited bubbles with smoothed surfaces cause a higher order phase transition. The order of 
the phase transition can be changed from a lower one to higher ones for the dilute and hot 
matter. 

The phase transition is understood in the following way. Hagedorn bubbles are suppressed 
strongly in the dilute hadronic matter because the pressure of the hadronic mass spectrum 
particles is higher than the Hagedorn bubble's internal pressure. These bubbles appear, 
however, due to the high thermal excitations. When their internal pressure exceeds the 
pressure of the gas of hadronic mass spectrum particles and other Hagedorn bubbles they 
start to expand at temperature just below the critical one. This expansion process depends 
on Hagedorn bubble's volume fluctuation and its internal color-flavor symmetry as well. The 
bubbles with lower internal color-flavor symmetries are easily excited than those with more 
complicated internal color-flavor structures. On the other hand, Hagedorn bubbles become 
more dominant for the large baryonic chemical potential. When the hadronic gas is heated, 
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the system passes a phase transition to quark-gluon plasma. 

Naively, the bubble's volume fluctuation depends on the reaction of the bubble's con- 
stituent quarks and gluons to the medium. It is expected the bubble's volume fluctuation 
becomes stronger in the dilute hadronic matter in particular near the critical temperature 
due to the high thermal excitations of the constituent particles. On the other hand, this 
volume fluctuation is supposed to be suppressed in the highly compressed matter in order 
to reduce the overlap effect among the hadrons. Two models for Hagedorn bubble's volume 
fluctuation have been considered for bubbles with various internal color-flavor symmetries. 
The first model is the Gaussian volume distribution function while the second one is the 
maximal volume distribution function, e.g. smearing the delta function 5(v — vq) to volume 
independent function. It is shown that the strength of volume fluctuation is essential for 
a hadronic phase consisting bubbles with low internal color-flavor symmetries to undergo a 
higher order phase transition while the complexity of the bubble's color-flavor symmetry is 
essential to determine the shape of the first order phase transition line for the large bary- 
onic chemical potential. The variation of the exponent 7(^5) in the medium for Hagedorn 
bubbles with specific internal color-flavor structures and the co-existence of the tri-critical 
point will be considered in the forthcoming work. 
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FIG. 1: The phase transition diagram in the baryonic density and temperature (ps — T) plane. The 
baryonic density is calculated for the hadronic phase below the phase transition line while for QGP 
is calculated above the phase transition line. The low thin lines show the phase diagram calculated 
using the Gibbs construction with various hadronic excluded volume. In Gibbs construction, 
Hagedorn bubbles (i.e. hadronic bubbles) are not included in the hadronic phase. The thick lines 
show pg G for the hadronic gas including Hagedorn bubbles with initial radii rpg = 1.0 fm. The 
Hagedorn bubbles' internal color-flavor structure is given by Eq.(172). The bag constant is taken to 
be i? 1 / 4 = 210 MeV. The above thin lines show p® GP for the QGP phase above the phase transition 
line for hadronic matter consisting Hagedorn bubbles with various internal color structure. 
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FIG. 2: Same as Fig. 2 but in the baryonic chemical potential and temperature — T) plane. At 
low temperatures, the hadronic matter consisting both the hadronic mass spectrum particles and 
Hagedorn bubbles shift the phase transition line to larger chemical potentials than hadronic matter 
consisting only the hadronic mass spectrum particles. The bag constant for Hagedorn bubbles is 
taken to be B 1 ^ = 210 MeV. The critical temperature T c (fj,B = 0) is found sensitive to B. 
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FIG. 3: The baryonic densities for both Hagedorn bubbles (i.e. hadronic bubbles) and the hadronic 

mass spectra particles in the hadronic phase below the phase transition line. The above thick 

oc p 

line shows the baryonic density for the QGP phase above the phase transition line for 

hadronic matter consisting hadronic mass spectrum particles and Hagedorn states with initial 
radius rpg = 1.0 fm for Hagedorn bubble's volume fluctuation. The excluded volume for the 
hadronic mass spectrum particles is taken to be rjj = 0.25 fm. At high temperatures and low 
chemical potentials, the hadronic mass spectrum particles is dominated while Hagedorn bubbles 
become the dominant for large chemical potentials and low temperatures, a) versus baryonic 
chemical potential [is- b) versus temperature. 
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FIG. 4: The phase transition diagram to QGP in the baryonic chemical potential and temperature 
(hb — T) plane for the hadronic phase consisting all the known hadronic mass spectrum particles as 
well as Hagedorn bubbles (i.e. hadronic bubbles) with various color-flavor internal symmetry and 
a specific volume fluctuation as in the model (I) given in the text. The phase transition diagram 
is displayed for various values of the color-flavor structures as determined by the pre-exponential 
power 7(a) factor 0.5, 1.00, 1.50, 2.00 and 3.00. The hadronic phase is chosen to be an ideal gas of 
the hadronic mass spectrum particles with the excluded volume of r# =0.30 fm as well as Hagedorn 
bubbles with an initial radii of rpg=1.0 fm. The density of states for Hagedorn bubbles is given 
by Eq.(172). The small excluded volume component for the hadronic mass spectrum particles and 
the large excluded volume component for the Hagedorn states approximation is considered in the 
numerical calculations. 
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FIG. 5: The pressure for the gas of hadronic mass spectrum particles and the pressure for the gas of 
Hagedorn states (i.e. hadronic bubbles) versus the baryonic chemical potential [Ib at temperatures 
just below the phase transition line. These pressures are displayed for Hagedorn bubbles with 
various values of internal color-flavor structure a = 47 — 1 factor of 0.5, 1.00, 1.50, 2.00 and 3.00 
where the density of states is given by Eq.(172). It is shown that the gas pressure for Hagedorn 
bubbles becomes larger than that for the hadronic mass spectrum particles for large baryonic 
chemical potentials. The small excluded volume component for the mass spectrum particles and 
the large excluded volume component for Hagedorn states approximation is considered in the 
numerical calculations. 
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FIG. 6: The sketch for the phase transition diagram outlining the shape and order of the phase 
transition between the hadronic and quark-gluon plasma phases. The hadronic gas consists all of 
the known hadronic mass spectrum particles and Hagedorn bubbles with various values for color- 
flavor correlations as determined by the pre-exponential power factor a = 47 — 1. The sketch is for 
model (I) with the density of states given by Eq.(172) (e.g. Hagedorn bubbles with strong volume 
fluctuation). The strength of the color-flavor correlation factor 7 = 7(^) (e.g. ct(fJ-)) changes the 
order of the phase transition for small baryonic chemical potentials and high temperatures while 
it modifies the shape of the phase transition at large chemical potentials and low temperatures. 
The phenomenological assumption is that the exponent 7 is modified self-consistently with respect 
to the baryonic chemical potential [ib- The variation of the phenomenological exponent 7 in the 
medium causes the existence of critical point in the phase transition diagram. 
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FIG. 7: Same as Fig. 4 but for density of states including the volume fluctuation consistently as 
given in model (II) (see Eq.(186)). 
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FIG. 8: The same as Fig. 5 but for the density Eq.(186), i.e. model (II). 



67 



Quark Gluon Plasma (Deconfined) 



3rd order 



2nd order (5/8>y>l/2) 
y=l/2 



1st order 



y>3/2 



Quark Gluon Droplet 

(Definite internal symmetry) 



Hadronic matter 




Y(M<b) 



B 



FIG. 9: Sketch of the phase transition diagram as given in Fig. 6 but for density of states for model 
(II) as given by Eq.(186) (e.g. soft volume fluctuation). In this scenario the color-flavor correlation 
is determined by the power factor a = 47 — ~. The points 7 = 3/2 and 1/2 correspond to a second 
and a third order phase transition, respectively. A first order phase transition takes place only for 
7 > 3/2. 
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FIG. 10: Rough Sketch of the dependence of Hagedorn bubble's phenomenological internal struc- 
ture parameter 7 = 7(/x_b) on the volume, the temperature (or the energy in the microcanonical 
ensemble) and the chemical potential. 
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